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FOREWORD 


This  report  was  prepared  by  John  W.  Lincoln  from  Structures  Division  of 
the  Directorate  of  Flight  Systems  Engineering,  Aeronautical  Systems 
Division,  Wright-Patterson  Air  Force  Base  Ohio.  Its  purpose  is  to  provide 
an  analytical  method  for  the  accurate  and  rapid  calculation  of  the  loads  on 
an  aircraft  during  landing.  The  report  is  written  in  two  volumes.  In  the 
first  volume  the  equations  of  motion  are  derived  and  in  the  second  volume 
the  computer  program  that  was  developed  from  these  equations  of  motion  is 
documented. 
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INTRODUCTION 


The  analysis  of  an  aircraft  making  a  landing  has  interested  dynamicists 
for  many  years.  This  interest  is  derived  from  the  technical  challenge  of 
the  problem  and  from  the  fact  that  the  landing  impact  is  an  important  loads 
source  for  both  land  and  carrier  based  aircraft.  The  inherent  nonlinear 
nature  of  the  equations  of  motion  requires  that  solutions  be  obtained  by 
numerical  integration.  The  challenge,  therefore,  is  in  establishing  the 
equations  such  that  solutions  may  be  obtained  with  speed  and  accuracy. 

Speed  is  required  because  of  the  need  to  minimize  the  computer  costs  derived 
from  the  large  number  of  cases  that  must  be  examined.  Accuracy  is  essential 
because  it  will  reduce  the  need  for  ground  testing  to  establish  the  forces 
and  the  energy  absorbing  capability  for  a  given  landing  gear  geometry. 
Further,  and  likely  more  important,  an  accurate  analysis  permits  one  to 
extrapolate  test  conditions  to  design  conditions  that  are  difficult  or  even 
hazardous  to  achieve  by  flight  test. 

The  landing  impact  analysis  in  this  report  is  believed  to  offer  both 
speed  and  accuracy.  Speed  is  obtained  through  taking  into  account  the 
compressibility  of  the  hydraulic  fluid  in  the  shock  strut.  The 
compressibility  is  accounted  for  in  the  equations  of  motion  through  the 
introduction  of  an  additional  first  order  differential  equation.  This,  in 
effect,  allows  the  integration  to  proceed  more  rapidly  than  that  obtainable 
if  the  ususl  incompressible  hydraulic  fluid  assumption  is  used. 

Accounting  for  the  compressibility  of  the  hydraulic  fluid  enables  the 
analyst  to  improve  the  accuracy  of  the  solution.  In  many  shock  struts, 
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particularly  in  articulated  landing  gears,  the  pressures  are  large  enough  to 
be  significantly  affected  by  the  hydraulic  fluid  compressibility.  The 
details  of  the  shock  strut  equations  that  account  for  compressibility  are 
found  in  Appendix  B.  Other  features  in  this  derivation  that  improve  the 
accuracy  of  the  solution  are  six  rigid  body  degrees  of  freedom,  flexible 
body  degrees  of  freedom  for  the  airframe  and  repr^oentation  of  the 
articulated  kinematics  of  the  landing  gears  in  combination  their  flexible 
degrees  of  freedom. 

The  equations  of  motion  are  derived  with  the  stroke  function  of  each  of 
the  shock  struts  used  as  a  generalized  coordinate.  This  feature  is 
especially  helpful  when  the  gear  is  articulated.  In  this  case  the  gear 
kinematics  and  gear  mass,  damping  and  stiffness  terms  become  functional 
dependent  on  the  shock  strut  stroke.  The  shock  strut  stroke  function  and 
its  derivative  are  obtained  directly  from  the  equations  of  motion  instead  of 
deriving  them  from  the  axle  displacemetit  and  velocity  functions. 

The  equations  of  motion  are  derived  from  Lag'^ange's  equations.  The 
rigid  body  motion  of  the  aircraft  is  expressed  in  terms  of  the  body  axis 
components  of  the  aircraft  velocity  vector  function  and  body  axis  components 
of  the  aircraft  angular  velocity  vector  function.  These  coordinates  are 
quasi-coordinates  and  the  usual  formulation  of  Lagrange's  equations  must  be 
modified  to  account  for  the  fact  that  these  coordinates  are  not  generalized. 
Lagrange's  equations  must  also  account  for  a  non-holonomic  constraint 
condition  on  the  vertical  motion  of  a  point  on  the  tire  footprint  in  contact 
with  the  ground.  This  constraint  condition  is  expressed  by  a  relationship 
involving  the  velocity  and  angular  velocity  terms.  It  is  found  that  the 
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Lagrangian  multiplier  in  the  equations  of  motion  is  the  vertical  ground 
reaction  on  the  landing  gear  tire.  The  modification  of  Lagrange's  equations 
fo.  quasi-coordinates  and  for  non-holonomic  constraints  is  given  in  Appendix 

C. 


The  aircraft  flexibility  is  represented  by  orthogonal  vibration  modes. 
The  influence  of  this  flexibility  is  included  in  the  calculation  of  the 
shock  strut  forces.  The  loading  on  the  aircraft  from  the  landing  impact  may 
be  calculated  from  the  displacement  or  acceleration  defined  by  the  flexible 
and/or  rigid  body  coordinates. 

Finally,  the  equations  are  included  for  the  calculation  of  the  forces 
from  a  cable  type  arresting  system.  It  is  assumed  that  the  aircraft  is 
equipped  with  an  arresting  hook.  The  equations  are  derived  so  that  the  case 
where  the  arresting  hook  meets  the  cable  off  center  and  the  aircraft  motion 
is  unsymmetrical  is  included.  The  condition  under  which  the  hook  will  or 
will  not  slip  on  the  cable  is  derived  and  the  resulting  cable  kinematics  are 
incorporated  in  the  equations.  The  details  of  this  derivation  are  given  in 
Appendix  F. 
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DERIVATION  OF  THE  EQUATIONS  OF  MOTION 


The  equations  of  motion  for  the  simulation  of  the  landing  impact  of  an 
aircraft  are  derived  below.  The  following  degrees  of  freedom  are  included 
for  the  airframe; 

Three  translational  degrees  of  freedom  for  the  rigid  airframe 
Three  rotational  degrees  of  freedom  for  the  rigid  airframe 
Flexible  body  degrees  of  f»'eedom  for  the  airframe  represented  by  its 
vibration  modes 

For  each  o*'  the  landing  gears  the  following  degrees  of  freedom  are 
included : 

Rigid  body  motion  of  the  gear  relative  to  the  airframe  that  is  dependent  on 
the  shock  strut  stroke  only 

Rigid  body  motion  of  the  gear  relative  to  the  airframe  that  is  dependent  on 
castoring  of  the  gear  (all  gears  are  assumed  to  have  this  degree  of  freedom) 
Flexible  body  degrees  of  freedom  for  the  gear  represented  by  its  vibration 
modes 

Rotational  motion  of  the  gear  wheel  about  the  gear  axle 

The  height  above  the  ground  of  the  ground  reference  point,  C  ,  in  the 

di 

wheel  midplane,  will  be  included  as  a  generalized  coordinate  in  the 
formulation  of  the  equations  of  motion.  However,  a  constraint  condition  on 
the  height  of  the  airframe  reference  point,  N,  above  the  ground  provides  the 
means  to  eliminate  this  coordinate  from  the  final  equations  of  motion. 
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The  six  rigid  body  motion  of  the  airframe  may  be  described  by  several 
different  coordinate  systems.  For  the  landing  impact  problem,  it  is 
convenient  to  use  the  body  axis  components  of  the  airframe  velocity  vector 
function  and  the  body  axis  components  of  the  airframe  angular  velocity 
vector  function.  In  general,  motion  associated  with  any  one  of  the  six 
degrees  of  freedom  can  have  an  influence  on  the  the  landing  gear  and  the 
airframe  loading. 

The  airframe  flexibility  may  have  a  significant  effect  on  the  shock 
strut  forces  and  consequently  must  be  included  in  the  equations  of  motion. 
The  airframe  flexible  body  equations  are  also  essential  if  the  airframe 
dynamic  response  loading  is  to  be  computed.  This  flexibility  may  be 
accounted  for  by  several  techniques,  howcer,  the  vibration  modes  which  are 
orthogonal  appear  to  be  the  most  efficient.  For  tha.s  derivation  it  is 
supposed  that  these  normal  modes  are  available.  It  is  further  supposed  that 
the  structural  damping  effects  can  be  simulated  adequately  by  linear  viscous 
damp-  ig  of  each  mod;,  independently. 

The  shock  strut  stroke  function  is  used  as  a  gencr.- zed  coordinate  in 
the  derivation.  Consequently,  ..he  shock  strut  force  functions  as  developed 
in  Appendix  B  are  expressed  in  terms  of  the  stroke  and  stroking  velocity 
functions.  Appendix  B  gives  the  neccessary  equations  for  two  of  the  many 
possible  gas  and  fluid  shock  strut  configurations.  For  both  of  these  shock 
struts  the  pressure  in  the  high  pressure  s'de  is  calculated  from  the 
compression  of  the  fluid.  At  a  time  t  the  rate  of  change  of  the  fluid 
compression  with  time  is  formulated  from  the  shock  strut  stroking  velocity 
and  the  voluae  of  fluid  passing  through  the  orifice.  Consequently,  another 


integr'^tion  is  added  to  calculate  the  compression  of  the  fluid. 

This  approach  may  appear  to  be  inefficient  when  compared  with  the  usual 
velocity  squared  damping  assumption.  However,  this  additional  integration 
provides  stability  to  the  numerical  integration  process  and  consequently 
permits  the  use  of  a  larger  step  size.  Treating  the  fluid  as  compressible 
also  provides  a  more  accurate  representation  of  the  shock  strut  fluid 
pressure  since  for  some  shock  struts  (particularly  in  articulated  gears)  the 
effect  of  compressibility  is  significant. 

There  are  several  unit  vector  systems  included  in  the  derivation  of  the 
equations  of  motion.  An  inertially  fixed  set  of  unit  vectors  is  used  as  a 
basis  from  which  the  motion  of  the  aircraft  is  defined.  At  a  time  t  after 
the  aircraft  has  contacted  the  ground  surface,  an  airframe  fixed  set  of  unit 
vectors  are  oriented  in  space  with  yaw,  pitch  and  roll  Euler  angles  for  the 
airframe.  To  provide  some  simplification  in  the  development  of  the 
geometric  data  and  mass  terms,  each  gear  has  a  set  of  unit  vectors  that 
permits  roll  and  pitch  of  the  gear  'eference  axes  relative  to  the  airframe. 
In  addition,  a  unit  vector  system  that  is  fixed  in  the  castored  gear 
reference  system  is  provided.  Since  the  wheel  plane,  because  of  geometrical 
constraints ,  may  not  be  convetiiently  oriented  by  the  gear  fixed  unit 
vectors,  a  separate  coordinate  system  is  fixed  in  the  wheel  plane.  These 
wheel-plane  unit  vectors  are  used  in  the  derivation  of  the  components  of  the 
ground  force  on  the  tire.  This  deri  atlon  is  given  in  Appendix  E. 
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Kinetic  Energy  Fo'^mulation 


The  derivation  of  the  equations  of  motion  is  accomplished  through 
Lagrange's  equations.  The  use  of  these  equations  is  believed  to  offer  a 
significant  advantage  over  the  direct  application  of  Newton's  second  law 
because  of  the  relative  ease  in  deriving  the  equations  for  the  articulated 
gear  geometries.  At  a  time  t  after  the  aircraft  has  contacted  the  groun^ 
surfac  the  kinetic  energy  of  the  aircraft  is  divided  between  the  airframe 
(i.e.  the  complete  aircraft  minus  the  gears)  and  the  gears.  The  gear 
kinetic  energy  is  divided  between  the  wheel  and  the  gear  structure  minus  the 
wheel.  The  following  definitions  are  needed: 

Q  is  an  inertially  fixed  ground  reference  point. 

N  is  a  reference  point  ir.  the  airframe.  The  derivation  is  made  with  the 
assumption  that  this  is  not  the  aircraft  center  of  gravity.  However,  if  the 
center  of  gravity  i^  used  for  N  then  some  simplifications  are  realized. 

T  is  the  a  gear  trunnion  to  airframe  attachment  point. 

Si 

A  is  the  point  that  is  common  to  ;.he  a  gear  wheel  raidplane  and  the  a  gear 
axle  centerline. 

Cp  is  the  point  that  ..s  contai.:ed  in  the  ground  surface  at  the  center  of 
a 

pressure  of  the  a  gear  tire. 

is  a  point  that  is  common  to  the  wheel  midplane  and  the  line  through  the 
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point  C  that  is  perpendicular  to  the  ground, 
r 

a 

H  is  the  point  of  attachment  of  the  arresting  hook  to  the  airframe.  It  is 
supposed  that  the  arresting  hook  rotates  about  an  axis  through  the  point  H 
and  parallel  to  the  lateral  axis  of  the  airframe. 

LP  is  the  lateral  pivot  point  in  the  arresting  gear  shank. 

P  is  the  aerodynamic  reference  point. 

TL  is  a  point  in  the  engine(3)  thrust  line. 

Vg  is  the  aircraft  volume  containing  a  set  of  mass  points  with  each  mass 
point  labeled  by  a  Cartesian  coordinate  system  fixed  in  the  airframe  which 
for  this  purpose  is  assumed  to  be  in  a  "jig"  or  undeformed  condition.  These 
mass  point  labels  stay  the  same  when  the  structure  is  deformed  as  a  result 
of  external  loads,  however,  in  this  case  the  mass  points  move  relative  to 
each  other.  The  concept  of  a  jig  condition  is  useful  for  the  definition  of 
certain  vectors  that  permit  the  rigid  body  motion  of  the  structure  to  be 
separated  from  the  motion  of  the  structure  due  to  deformation. 

N  is  the  number  of  gears  on  the  aircraft,  a  is  in  [1,  N„]  and  V„  is  the  a 
u  vj  G 

a 

gear  volume  excluding  the  wheel.  The  gear  mass  points  are  labeled  in  a 
manner  similar  to  that  used  in  the  airframe.  The  jig  condition  of  the  gear 
is  defined  with  the  shock  strut  in  the  fully  extended  position. 

is  the  a  gear  wheel  volume, 
a 
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Pn  is  the  sicple  aariace  such  that  if  (x„\  Xr,^>  is  in  V_, 

O  O  D  D  B 

P,(x„\  x„^,  x„^)  is  the  mass  density  at  the  airframe  point  labeled 

O  D  D  D 

,  1  2  3x 

(Xb  ,  Xb  ,  Xg  ). 


P  is  the  simple  surface  such  that  if  (x_  \  x„  x_  is  in  V_  , 
u  vi  G  G  G 

a  a  a  a  a 

12  3 

P,  (Xq  ,  Xq  ,  Xq  is  the  mass  density  at  the  a  gear  (without  wheel) 
a  a  a  a 

12  3 

point  labeled  (x_  ,  x„  ,  x_  ) . 

G  G  G 

a  a  a 

P  is  the  simple  surface  such  that  if  (Xy  \  ’ 

a  a  a  a  "a 

(Xy  \  mass  density  at  the  a  gear  wheel  point  labeled 

a  a  a  a 

,  1  2  3v 

*w  ’  * 

a  a  a 

-  12  3 

R„  is  the  vector  function  such  that  if  (x„  ,  x_  ,  x„  )  is  in  V_  and  t  >  0, 

D  D  D  D  D  — 

“12  3 

R„(Xn  ,  x_  ,  x„  ,  t)  is  the  vector  at  the  time  t  from  the  ground  reference 

r  D  D  D 

12  3 

point  Q  to  the  airframe  point  labeled  (x^  ,  x_  ,  x  ) . 

D  D  B 


-  12  3 

R^  is  the  vector  function  such  that  if  (x„  ,  x_  ,  x_  )  is  in  and 

u  u  Cj  G  G 

a  a  a  a  a 

12  3 

t  0,  R(,  (x^  ,  X(,  ,  Xq  ,  t)  is  the  vector  at  the  time  t  from  the  point 

a  a  a  a 

12  3 

Q  to  the  gear  (without  wheel)  point  labeled  (x^  ,  x_  ,  x^  ) . 

G  G  G 

a  a  a 


•  12  3 

Ry  is  the  vector  function  such  that  if  (Xy  »  Xy  ,  Xy  )  is  in  Vy  and 

a  a  a  a  "a 

-  12  3 

t  0,  Ry  (Xy  ,  Xy  ,  Xy  ,  t)  is  thc  vcctor  at  the  time  t  from  the  point 
a  a  a  "a 

12  3 

Q  to  an  the  gear  wheel  point  labeled  (x,,  ,  x,,  ,  x,,  ). 

WWW 
a  a  a 


9 


Therefore,  if  T  is  the  simple  graph  such  that  if  t  0,  T(t)  is  the 
kinetic  energy  at  the  time  t  of  the  collection  of  mass  points,  then  with  the 
derivative  notation  described  in  Appendix  A,  T  may  be  expressed  by 


The  Ro  vector  function  must  be  expressed  in  terms  of  constituent  vector 

D 

functions  in  order  to  develop  the  equations  of  motion.  For  this  purpose  the 
following  definitions  are  needed: 

r  is  the  vector  function  such  that  if  t  0,  r(t)  is  the  vector  at  the  time 

t  from  the  point  Q  to  the  jig  condition  location  of  the  airframe  reference 

1  2  3 

point  N  labeled  (x„  ,  x„  ,  x_  ). 

12  3 

L„  is  the  vector  function  such  that  if  (x^  ,  x_  ,  x^  )  is  in  V„  and  t  >  0, 

D  D  D  D  D  — 

-12  3 

L„(x„  ,  x„  ,  Xn  ,t)  is  the  vector  at  the  time  t  from  the  jig  condition 

D  D  D  D 

location  of  the  airframe  reference  point  N  to  the  jig  condition  location  of 

1  2  3 

the  airframe  point  labeled  (x„  ,  x„  ,  x_  ) . 

D  o  B 

12  3 

Uj,  is  the  vector  function  such  that  if  (x_  ,  x  ,  x^  )  is  in  V_  and  t  >  0, 
-12  3 

Ug(Xg  ,  Xg  ,  Xg  ,  t)  is  the  vector  at  the  time  t  from  the  jig  condition 

1  2  3 

location  of  the  airframe  point  labeled  (x  ,  x  ,  x„  )  to  the  actual 

B  B  B 

location  of  this  point. 
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Therefore,  the  vector 


V^b’-  *B^'  ''  =?(')■>  *B^’ 

*  "b'-'b'’  -b^’  ’'e^- 

can  be  used  to  derive  the  vector  function  R  . 

D 

The  vector  functions  for  expansion  of  are  defined  as  follows: 

U 

a 

is  the  vector  function  such  that  if  t  0,  (t)  is  the  vector  at  the 

a  a 

time  t  from  the  Jig  condition  location  of  the  airframe  reference  point  N  to 

the  jig  condition  location  of  the  point  T  . 

a 

-  12  3 

G  is  the  vector  function  such  that  if  (x-  ,  x„  ,  x„  )  is  in  V_  and 

a  u  G  G  u 

a  a  a  a 

-  1  2  3 

t  >  0,  G  (x-  ,  X-  ,  x-  ,  t)  is  the  vector  at  the  time  t  from  the  jig 

—  a  G  G  G 

a  a  a 

condition  location  of  the  point  T  to  the  jig  condition  location  of  the  a 

a 

1  2  3 

gear  point  (x„  ,  x„  ,  x„  ).  The  jig  condition  for  the  a  gear 

u  u  u 

a  a  a 

is  determined  with  the  shock  strut  stroke  equal  to  zero  (fully  extended 
gear) . 


s  is  the  simple  graph  such  that  if  t  >  0,  s  (t)  is  the  a  gear  stroke  at  the 
a  a 

time  t. 


— .  12  3 

S  is  the  vector  function  such  that  if  (x_  ,  x„  ,  x„  )  is  in  V„ 

a  G  G  G  G 

a  a  a  ; 

1 


and  t  >  0,  “  (x„  ' ,  x_  x  s  ( t) ,  t)  is  that  part  of  the  vector 
—  a  u  vj  u  a 

a  a  a 


at  the  time  t  from  the  jig  condition  location  of  the  a  gear  point 


12  3 

(Xq  ,  Xj,  ,  Xq  )  to  the  actual  location  of  this  point  which  is  derived 
a  a  a 

from  the  rigid  body  motion  of  the  gear  as  result  of  the  stroking  of  the 
shock  strut. 


is  the  simple  graph  such  that  if  t  ^  0,  0^  (t)  is  the  a  gear  castor 
a  a 

angle  at  the  time  t.  Normally,  only  th  auxiliary  gear  is  permitted  to 

castor.  However,  for  symmetry  of  the  equations  of  motion,  all  of  the  gears 


in  this  derivation  will  be  assumed  to  have  this  degree  of  freedom. 


-  12  3 

U  _  is  the  vector  function  such  that  if  (x„  ,  x_  ,  x„  )  is  in  V-  and 

DU  G  G  G  G 

a  a  a  a  a 

-  12  3  /3 

t  2.  0,  UgQ  (Xq  ,  Xq  ,  Xq  ,  a^(t),  C/q  (t),  t)  is  that  part  of  the  vector 
a  a  a  a  a 

at  the  time  t  from  the  jig  condition  location  of  the  a  gear  point 
1  2  3 

(Xq  ,  Xq  ,  Xq  )  to  the  actual  location  of  this  point  which  is 
a  a  a 

derived  from  the  deformation  of  the  airframe.  Note  that  the  shock  strut 
stroke  and  the  castor  angle  dependence  appears  explicitly  in  this  vector. 
This  provides  for  a  simpler  formulation  when  the  deformation  vector  is 
expressed  in  terms  of  its  components. 


-  12  3 

U„  is  the  vector  function  such  that  if  (x„  ,  x„  ,  x„  )  is  in  V_  and 

u^  G  G  G  G 

a  a  a  a  a 

t  0,  Uq  (Xq  \  *G  that  part  of  the  vector 

a  a  a  a  ^a 

at  the  time  t  from  the  jig  condition  location  of  the  a  gear  point 
1  2  3 

(Xq  >  *G  ’  *G  ^  actual  location  of  this  point  which  is 

a  a  ^a 

derived  from  the  deformation  of  the  a  gear. 


Therefore,  the  vector 


(3) 


(x-  \  x_  X  t)  =  r(t)  +  i„  (t) 

u  U  vj  i 

a  a  a  a  a 

-  1  2  3  <—123 

+  G  (x„  ,  x_  ,  x„  ,  t)  +  r;  (x„  ,  x_  ,  x_  ,  3  (t),  t) 

a  u  u.  1.  I-  I.  o 

a  A  a 


'a  G  ’  G  ’  G  ’  a 
a  a  a 


^BG  ^*G  ’  *G  ’  ’'g  ^G 

a  a  a  a  a 

+  Uq  (Xq  ,  Xq  ,  Xq  s^(t),  (t),  t) 

a  a  a  a  a 


can  be  used  to  derive  the  vector  function  R„  . 

u 

a 


The  final  vector  function  needed  for  the  kinetic  energy  formulation  is 

R^  .  The  vector  functions  used  to  define  R^  are  defined  as  follows: 

”a  a 

1  2  3 

If  (x_  ,  x_  I  X-  )  is  the  label  of  the  a  gear  axle  point,  A  ,  and 

°A  °A  °A  ^ 

a  a  a 

if  t  2.  0,  then 

-  -  1  2  3 

8»  (t)  =  Gg^(xQ  ,  Xq  ,  Xq  t), 

“a  AAA 

a  a  a 


(s  (t),  t)  =  s  ( Xp  ,  x_  ,  x_  ,  s  (t),  t) , 


'A  '  a' 
a 


’a'^'G,  ’  '‘G,  ’  G,  ’  a' 

AAA 
a  a  a 


^BG.  ^G  ■  ^BG  ^’'g  ’  *G  ’  *G  ^G  ’ 

A  a  a  A  A  A  a 

a  a  a  a 


U  (s  (t),  ^  (t),  t)  =  U  (x  \  X  X  3  (t),  (t),  t), 


a  A  A  A 

a  a  a 


at  the  time  t. 


-  12  3 

W  is  the  vector  function  such  that  if  (x„  ,  x  ,  x  )  is  In  V  and 

3  w  w  w  w 

a  a  a  a 

-  1  2  3 

t  2  G»  ’  ’'w  ’  ®  vector  at  the  time  t  from  the 

a  "a  ”a 

axle  point  A  to  the  jig  condition  location  of  the  point  in  the  a  gear 

3 


labeled  (x^  \  x^^  x^  ^). 

a  a  a 
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Therefore,  the  vector 


a  a  a  a 


+  I_  (t)  +  g.  (t)  +  (s^(t), 


A  '  a' 
a 


t)  (4) 


+  0  (s  (t),  0  (t),  t)  +  U_  (s  (t),  Q  (t),  t) 


BG.  'a 
A 

a 


"a  ■  ^ 

a 


+  W^(Xy  \  Xy  3^(t),  t) 

a  a  a 


can  be  used  to  derive  the  vector  function 


To  be  useful  for  the  kinetic  energy  equation  (equation  (1)),  the 
structural  deformation  vectors  appearing  in  equations  (2)  through  (4)  must 
be  expressed  in  terms  of  generalized  coordinates.  This  may  be  done  through 
the  use  of  the  vibration  modes  of  the  structure.  Orthogonality  of  these 
modes  is  assumed  in  this  derivation  since  this  property  is  typically  used  in 
practice. 


Suppose  that  N  airframe  vibration  modes  are  used  to  define  the 

otL 

airframe  deformation.  Further  suppose  that  b  is  an  integer  in  [1,  N__]  and 

BE 

that  q  ^  is  the  simple  graph  such  that  if  t  >  0,  q_*^(t)  is  the  bth  airframe 

D  —  D 

vibration  mode  displacement  at  the  time  t. 

Also,  suppose  that 

<J)_  is  the  vector  function  such  that  if  (x„  ,  x_  ,  x„^)  in  V„  and  t  >  0, 

£3^  D  B  B  B  — 

T  f  V  ^  2  3 

B  »  ‘  airframe  vibration  modal  vector  at  the 

12  3 

time  t  for  the  airframe  point  labeled  (x„  ,  x„  ,  x„  ) , 

B  B  B 
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""  12  3 

Ct)„_  is  the  vector  function  such  that  if  <x-  ,  x„  ,  X-  is  in  V. 

^  dU  u  U  vj  L 

a  a  a 


t  >_  0,  (Xq  \  *G 

a.  a  a  a  a 

D 


the  bth 


airframe  vibration  modal  vector  at  the  time  t  for  the  a  gear  point  labeled 
/  ■>  2  3. 

(Xg  ,  Xq  ,  Xq  ), 

a  a  a 


<i>  „„  is  the  vector  function  such  that  if  t  >  0, 


^BG  (s^(t),  0^  (t),  t)  is  the  bth  airfr. 

a 

^b 

at  the  time  t  for  the  a  gear  axle  point  A  . 

a 


the  bth  airframe  vibration  modal  vector 


Now  suppose  that  for  the  a  gear  there  are  N-„  gear  vibration  modes 

uo 

a 

used  to  define  its  deformation.  Suppose  that  d  is  an  integer  in  [ 1 ,  ] 

d  d 

and  that  q_  is  a  simple  graph  such  that  if  t  >  0,  q_  ^(t)  is  the  dth  a 

li  —  u 

a  a 

gear  vibration  mode  displacement  at  the  time  t.  Also  suppose  that 


12  3 

_  is  the  vector  function  such  that  if  (x_  ,  x  ,  x„  )  is  in  V  and 


G  ’  G  ’  G 
a  a  a 


12  3  ZD 

t  0,(J)q  (Xq  >  *  *G  ^G  ^  gear  vibration 

a.  a  a  ^a  a 

d 

12  3 

modal  vector  at  the  time  t  for  the  a  gear  point  (x„  ,  x„  ,  x_  ), 

u  u  0 

a  a  a 


Ct)-  is  the  vector  function  such  that  if 
H^Ga 

®d 


t>0,  (sjt).  Oq  (t),  t) 

^A  ^a 

^d 


is  the  dth  a  gear  vibration  modal  vector  at  the  time  t  for  the  a  gear  axle 


point  A^. 
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Therefore,  the  deforoation  vectors  may  be  expressed  as  follows: 

r,  ^  1  2  3  X  /  1  2  3  b... 

^  "^B  *B  ’  *B  ’  ’^B  ’  ^  ’ 

b 

where  summation  on  the  index  b  is  implied. 

^BG  ’  *G  ’  ’'g  ^G  " 

a  a  a  a  a 

<i^BG  ^"'g  ""g  ’^G 

a  a  a  a 

b 

OfiG  (s^(t),  0^  (t),  t)  =(i)gQ  (s^(t),  0Q  (t),  t)  qg^(t). 
a  A  a 

^  % 

®G  ^*G  ’  *G  ’  *G  ^G  ' 

a  a  a  a  a 

i^G  ^"'g  *G  *G  ^G  ‘^G 

a.  a  a  a  a  a 

a 

(3^(t),  0^  (t),  t)  =^Q  (3^(t),  0q  (t),  t)  Qq  ‘^(t). 

’a  “  a  A„  a  a 

®  d 

At  this  point  in  the  derivation  consideration  must  be  given  to  the 
order  of  magnitude  of  the  deformation  of  the  structure.  The  general  case  of 
large  deformations,  which  is  described  in  Reference  (1),  could  be  applied. 
However,  this  additional  complication  for  most  aircraft  structures  does  not 
appear  justified.  Therefore,  it  will  be  assumed  that  the  deformation  is 
small  enough  such  that  terms  in  the  kinetic  energy  expression  involving  the 
functions  and  the  q_  functions  can  be  eliminated.  Also,  it  will  be 

D  U 

a 

assumed  that  the  motion  of  the  body  axis  reference  axes  is  completely 
defined  by  the  rigid  body  equations  of  motion.  This  assumption  is  also 
justified  from  the  assumption  of  small  structural  deformations. 

Therefore,  from  equations  (2),  (3)  and  (4)  the  velocity  vectors  may  be 
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written  with  the  aid  of  the  following  definitions: 


is  the  vector  function  such  that  if  t  >  0,  is  the  airframe 

O  —  D 

angular  velocity  vector  at  the  time  t. 


is  the  vector  function  such  that  if  t  >  0,  (t)  is  the  a  gear 

u  —  u 

a  a 

castoring  angular  velocity  vector  at  the  time  t. 


(t)  is  the  vector  function  such  that  if  t  0,  (t)  is  the  a  gear 

a  a 

wheel  angular  velocity  vector  at  the  time  t. 


Thus,  with  the  assumption  of  small  deformations,  the  velocities  may  be 
written  in  the  following  vector  format: 


*B^’  *B^’  +i'^g(t)  X  Lg(Xg\  x^^,  x^^,  t) 

t 

+  C^Q  (Xg  ,  Xg  ,  Xg  ,  t)  Qg  '(t). 

b 


Rq  .j  (Xq  \  Xq  Xj,  t)  =  r'(t)  +^lg(t)  X  (t) 
a  ’  t  a  a  ^a  a 

+  (^g(t)  +  Oq  (t))  X  \  Xq  Xq  t) 

a  a  a  a 

(t))  x  ~  (x„  \  X  X  s  (t),  t) 
o  j  a  u  u  u  a 

a  a  a  a 

^a;a  ’^G 

a  a  a  a 

■^^BG  ^*G  ’  *G  ’  *G  ’  ^G  '^B 

a.  a  a  a  a 

D 

^G  ^*G  *G  ’'g  ‘^G 

a .  a  a  a  a  a 

d 
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a  t  a  a  a 


t)  =  ?’(t)  X  (t) 


+  (i7_(t)  +  Q,„  (t))  X  g  (■ '  +  (i7p(t)  +Q.  (t))  X  'b  .5  (t),  t) 

DU  A  DU  rt  a 

_  a  a  a  a 

+  (3  (t),  t)  s  '(t) 


’A  'a' 
a;s 

a 


BG,  a’ 
A 


+  <J)p  (3  (t),  Cn 


'•'a''"'  '^G  ^G 

A  a  a 

+  (i7g(t)  +  n.Q  (t)  (t))  X  W^(Xy  \  x^^  s^(t),  t) 

a  a  a  "a  a 

^  \;3  *W 

a  a  a  a 


The  vector  functions  Rq,  ,  R^,  ,j  and  R  may  be  easily  derived 

“’•‘’t  a’^t  "a’'^t 

from  equations  (5),  (6),  and  (7). 


Before  Rg.j  >  Bq  and  ^  .j.  are  substituted  int-o  the  kinetic  energy 
’■^t  a’^^t  a’'^t 

equation  (equation  (1))  it  is  convenient  to  express  the  vectors  appearing  in 
these  functions  in  terms  of  their  components.  To  accomplish  this,  several 
sets  of  right  handed  orthogonal  unit  vector  functions  need  to  be  defined. 

Suppose  each  of  !„  ;  b  =  1,  3  is  a  unit  vector  function  such  that  if 
®b 

t  >  0,  (t)  is  the  vector  parallel  to  the  airframe  longitudinal  axis  and 

directed  forward,  I  (t)  is  parallel  to  the  airframe  lateral  axis  and 

directed  to  the  right  relative  to  the  pilot  and  I  (t)  is  parallel  to  the 

airframe  vertical  axis  and  directed  down  relative  to  the  pilot  at  the  time  t. 


To  orient  the  airframe  in  space  relative  to  an  inertial  reference  frame 


defined  by  the  unit  vectors  b  =  1 
For  this  purpose  suppose  that  each  of 


3,  a  set  of  Euler  angles  is  used, 
and  is  a  simple  graph  such 
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that  if  t  0,  4'm.  0m  and  are  at  the  time  t  the  yaw,  pitch  and 

roll  tuler  angles  of  the  airframe.  These  Euler  angles  are  used  in  the 
transformation 


Ib  I.,’ 

C  C 

where  the  components  of  are  derived  in  Appendix  G. 

D 

C 

From  Appendix  G  it  is  seen  that  and  are  in  the  ground  plane  and 
is  perpendicular  to  the  ground. 


Also,  suppose  that  each  of  i  ;  b  =  1,  3  is  a  unit  vector  function 

U 

^  r 

and  each  of  ij  ^  and  is  a  number  such  that  7^  ^  is  the  pitch  Euler 

a  a  '^a 

angle  and  is  the  roll  Euler  angle  for  the  a  gear  relative  to  the 

airframe  such  tnat  if  t  0,  each  of  i^  (t);  b  =  1,  3  are  fixed  relative 


to  the  a  gear  at  the  time  t  and  are  oriented  relative  to  the  airframe  fixed 
unit  vectors  through  the  transformation 


^G 

a  a  b 

c  c 

whore  the  components  of  ^  are  derived  in  Appendix  G  in  terms  of  the 

Euler  angles  Tj ^  and  . 

a  a 


Suppose  that  each  of  ;  b  =  1 ,  3  is  a  unit  vector  function  such  that 

% 

if  t  0,  each  of  (t);  b  =  1,  3  is  a  vector  fixed  in  the  a  gear  at  the 
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tiue  t. 


^3  used  to  orient  these  vectors  with  respect 


The  castor  angle,  0 


to  i„  (t);  b  =  1,  3  vectors  through  the  transformation 
G 

% 

Jg  =  Qg  ‘g  Sbe  S°*' 

a  a 

b  g  c 

A  g 

where  Q„  (t)  is  defined  in  Appendix  G. 
a 

g 

It  is  also  useful  to  define  the  transformation  from  the  airframe  fixed 

unit  vectors  I  (t);  b  =  1,  3  to  the  a  gear  fixed  unit  vectors 
®b 

(t);  b  =  1,  3.  This  may  be  done  through  a  combination  of  the 
G 

% 

transformations  defined  above  as  follows: 


A  e ,  ^ ,  o  /  d 


(t)  S„„  5^^ 


Jg  (t)  =a„  ’(t)  “(t>  i  (t)  0,, 

a.  a  a  d 

b  g  c 


=  a,  "(t)  i  (t). 


The  components  of  Q  a>’e  given  in  Appendix  G. 

G 


Suppose  that  each  of  ;  b  =  T ,  3  is  a  unit  vector  function  such  that 

U 

^b 

if  t  >  0,  each  of  k  (t);  b  =  1,  3  is  fixed  in  the  a  gear  wheel  at  the  time 
—  G 

^b 

t  and  k_  (t)  is  a  unit  vector  parallel  to  the  a  gear  wheel  axle. 

G 

^2 

Therefore,  the  vector  k„  (t)  may  be  expressed  in  terms  of  the  vectors 

G 

^2 

J_  (t);  b  =  1,  3  by  the  relation 
G 

a 


i^G  ( t /  =  /^G  ^ ^ )  Jg  ’ 

^2  ^2 

where  from  the  gear  kinematics  there  exists  a  set  of  functions 


;  b  =  1 ,  3  such  that 

U 

^2 


/3q  ^(t)  =  Bg  ‘'(s^(t));  b  =  1,  3. 


Now  suppose  that  each  of*  ly  ;  b  =  1 ,  3  is  a  unit  vector  function  such 

®b 

that  if  t  0,  (t)  is  (t)  is  orthogonal  to  I^g  (t)  and  (t), 

^3  _  _  ^2  ^3 

and  ly  (t)  is  orthogonal  to  (t)  and  (t)  at  the  time  t. 

^2  ^3 

The  transformation  from  b  =  1,  3  to  (t);  b  =1,  3  is  defined  in 
Appendix  E  and  is  given  by 


iu  (t)  =  X  I  5.  8°®, 

W  ^  W  c  be 


b  g 

With  the  aid  of  the  unit  vector  functions  defined  above  the  linear  and 
angular  velocity  functions  are  defined  by 


r  ’  -  V  ^  I 

®  X’ 


^B  ^  ^B  ^B  ’ 
0 


'^G  ^G  ’ 
a  a  a^ 


n y  - 

a  a  a. 
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The  vector  functions  in  equations  (5),  (6)  and  (7)  may  be  defined  by 


Finally,  the  modal  vector  functions  for  the  deformation  of  the  airframe 
and  the  gears  are  defined  by 


^  BG  ■  BG 

b  be 


BG 


■  ^BG  ^G 
A  A  a 

S  ° 


^G  ■‘i^G  ^G  ’ 

b  be 
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Therefore,  with  these  definitions  the  kinetic  energy  function  T  may  be 


expressed  in  the  following  form: 


Vg  Vg  Hg  Hg 

I  =  <»  b  o>  C  V  »  <«  b 


)[s^,  ;  a  =  1,  (8) 


Nq  n,Q 


+  1/2  5;  (M  ^  ^)[s J  (^2^  +  1/2  2  (M 

a=1  ^  °a  a=1 


Nq 

2  .  1/2  I  (M  ^  ^  (H, 

a=1  a 


N„  3  3  c 

G  a  a 

+  1/2  T  (M  )[3  ]  (s  ' )  +  1/2  (M 

a  a 
a=1 


‘’b 


)[s^,  ;  a  =  1,  Ng]  qg^'  qg^"* 


b  c'^  a’ 


.  1/2  i  (M  W)[s^,  ]  q^"-  q^"* 

a=1  a  a  a 


B  '^‘'B 


)t3^,  ;  a  =  1,  Nq]  Vg  Og 


b  c'"  a’  ^G 


'^B  ^G 


Ng  Vg  s^ 

^  h  ^G  ^  ^G  ^  b  ^^®a’  ^G  ^  '^B  ®a' 

a  a  as1  a 


^B  ^^B 


'^G  '^B  <’g 


b  o^^®a’  ^G  ’  ®  ^  ^  b  c^'-^a’  ^G  ^  ^R  *^G  ' 


•G"  B  '  b  c'^  a’  G  "  B  '’G 

a=1  a  a 


Ng  rigH^ 


*'*G  ^  B  ^  W 

i  («  'b  '")[=,.  @0  *  2  <H  b  "ns^.  0^:a^'‘n„ 

a  a  a=  1  a  a 


Nf,  Hr  s  Hr  'Ir 

I  («  ,  )ts^,  0^  in/  ♦  (M  ^ 

a=  1  a 


)[s^,  0r.  >  ^ 


b  c ' "  a  ’  G 


2  <"  b  \»=b’  00  ine”  c,' 

a=  1  a  a 


^G  U  ^G  G  ®a 

+  Z  (M  ^  ^)[s^]f7g  n  „  +  2  (M  ^  )[s3]nc  s  ' 

a=1  a  a  a=1  a 
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Nq  Qg 

-  X  (M  ^ 
a=1 


Nq  Hq  Qq 

^  *  '*’  ^  b^'-^a’  ^G  ^  ^G  ‘^G 

a  a  a=l  a  a  a 


"g  =a  ‘’■B 

Z  («  I 

a=l 


"o  Ig 


b>t=a-  &C  '  =a'  "b'’'  *  S 

a  a=1 


b>‘\'  ^G  J  =a’  '>0  '’' 
a  a 


"g  ’b  ‘>g 
Z  <«  b  ' 

a=1 


Or,  ^  ^r, 


b  c'"  a’  '-'G  "  '^B  ^G 
a  a 


where  the  mass  terms  in  the  expression  above  are  defined  in  Appendix  D  and  the 
bracket  product  notation  is  defined  in  Appendix  A. 


Suppose  that  K  is  a  simple  surface  such  that  the  kinetic  energy  function  T 
given  by  equation  (8)  may  be  written  in  the  form 

T  =  K[Vg\  Vg^  Vg^  Hg^,  Hg^,  Hq  ,  ,  Hq  ,  (9) 

1  NG  ^ 

n  n  s'  s  '  a  Q  Q 

Oq  f  ••  •  »  Oq  ^  ••  •  . 


It  is  seen  that  the  function  T  formulated  in  equation  (9)  is  expressed 
partially  in  terms  of  coordinates  (i.e.  functions)  referred  to  in  Reference  (2) 
as  quasi-coordinates.  Namely,  the  body  axis  components  of  the  airframe  velocity 
vector  function  and  the  body  axis  components  of  the  airframe  angular  velocity 
vector  function  are  quasi-coordinates.  Appendix  C  describes  the  modifications 
to  Lagrange's  equations  when  quasi-coordinates  are  used. 


Suppose  that 


.  =  1, 


3  are  the  components  of  r  (the  vector  from  the 
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point  Q  to  the  jig  condition  location  of  the  point  N)  in  the  inertial  reference 
frame  defined  by  b  =  1,  3-  It  is  useful  (e.g.  in  Appendix  F)  to  use  the 
symbols  d,  s  and  v  as  alternate  designations  for  these  components. 


Therefore,  if  'P.  9  and  are  the  yaw,  pitch  and  roll  Euler  angle 
functions  which  are  discussed  in  Appendix  G,  then  consistent  with  the  notation 
of  Appendix  C,  the  following  generalized  coordinates  are  used  for  the  airframe 


rigid  body  motion: 

2  r  2 

q  =  L  =  s. 

=  V. 


-.4,. 
Q®  -  4. 


Also,  the  following  relations  are  consistent  with  Appendix  C  notation; 


-  V  ^ 

■  ''b  • 

2  2 

=  V32. 

_  V  3 
V  -  Vb  . 


v^-Hb'- 


The  transformation  from  the  quasi-coordinates  to  the  generalized 
coordinates  is 


where  from  Appendix  G  it  is  found  that 


/S]  =  cos  [0]  co3[^],  -  -  sin  [^]  cos[C^]  +  =  sin[^]  3in[(p>]  + 

3in[(^]  sin[^]  cos  cosicj)]  sin[^]  cos  ['#'], 

=  cosl0]  3in[^3,  =  cos[(^]  cos  I'P:  *  -  -  sin[c/)]  cos  i4i  * 

3in[<^]  3in[^]  3in[T],  oo3[C^]  sin[^]  3in[^], 


25 


:  -  sin[^], 

=  0, 

Bl 

•  H 

=  1, 

=  0, 

and 

all  other  ^ 

^  =  0. 
b 

The  elements  of  the 

vS 

Qg  q^'  are 

a,' 

=  cos[^]  cos 

>c^]. 

=  -  sinl^l  cos[(^]  + 

sin 

i[(^]  sint0] 

cos[^] , 

a? 

=  sin[^]  sin[<^]  + 

co3[(j>]  sin[^] 

cos[^] , 

=  0, 

=  cosicf)], 

6 

14 

=  -  siniCp], 

and  all  other  Q.*  =  0. 


=  sin[(f>]  coslQ'i, 
=  oo3[Cj>-i, 

B\  =  tan[@]  sin[(^] , 
=  ain[<^]/cos[6^], 

inverse  transformation 


=  coa[^]  sin 

l^h 

CM  CM 

d 

=  oo3[(^]  cos 

sin 

i[<^] 

sin[0] 

sincV^] 

3in[(t]  cos[^] 

co3[<^]  sin[^] 

sin[^] 

a  5 

=  1, 

=  0, 

a« 

=  0, 

=  oosLcf^]  cos  [0], 

0l  =  -  sin[(^], 

-  tan[^]  co3[(p>], 
=  cos[<^]/cos[^] , 

=  -  sin[^], 

CL^  =  3in[<jb]  cos[^], 

CL^  =  C03[(j)]  ooslQ], 

Clg  =  -  sin[^], 

^6  =  cos  [0]  sin[cf>h 
=  co3[0]  cos  [</)], 
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For  a  in  [ 1 ,  6]  and 

b  in  [1,  6],  the  functions 

D 

which  are  defined 

i  of  a®  and  yg®  in 

equation  (C-11)  of  Appendix  C  are 

Q]  =  0, 

nj  =  0, 

n;  =  -  V33. 

-  Vb'. 

0.1  =  0. 

ol  =  -O^. 

-  V  ^ 

-  Vg  , 

ol  =  0. 

ol  --s'. 

0l--0^\ 

ol  -  0, 

03  -  _  Y  ^ 

- B  ’ 

Ol  =  V3'. 

ol  -  0. 

II 

o 

=  -Os\ 

o-l -a^\ 

ol  -  0, 

0^  - 0^  1 

ol  ^0^'. 

Of,  -  °> 

and  all  other  ^7^  =  0* 
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The  Forces  Acting  on  the  System 


The  quasi-coordinate  force  terms  may  be  defined  through  equation  (C-l8) 
of  Appendix  C.  Included  in  these  terms  are  the  body  gravitational  force, 
the  arresting  gear  force,  the  ground  force  (parallel  to  the  ground  plane), 
the  aerodynamic  force  and  the  thrust  force.  There  are  also  some  generalized 
coordinate  force  terms  that  must  be  included  in  the  equations  of  motion. 
These  include  the  shock  strut  forces,  the  airframe  modal  elastic  and  damping 
forces,  the  gear  modal  elastic  and  damping  forces  and  the  gear  stiffness  and 
damping  forces  associated  with  castoring.  The  notation  on  the  left  side  of 
equation  (C-l8)  will  be  used  for  both  quasi  and  generalized  coordinates. 


Consistent  with  the  definitions  given  in  Appendix  C  and  with  g  equal  to 
the  gravitational  acceleration  the  virtual  work  done  by  the  gravity  force  is 


a  a  a  a 


Therefore , 


Ss 


and 


based  on  Appendix  C,  the 
S Ry  may  be  replaced  by 


a  a 

Rp  and  R  respectively, 

a’  t  "a’^t 


vector 


functions 
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Thus, 


,(t)  = 


”xg''  ■  l3>  « 


ij\  -V  \Bd,, 


*  "w  ;i,  ’•  V  V  ^  Pm  b  ■“''w  ■ 

a=1<^W  a’tW  W  W  a  a 

a  a  a  a 


If  equations  (5),  (6)  and  (7)  are  substituted  into  equation  (10)  and  if 
the  mass  terms  defined  in  Appendix  D  are  utilized,  then  the  following 
functions  are  found  to  be  the  coefficients  of  the  coordinate  velocities: 


"b  'b  *6  ^  ^  , 

S„  (t)  =  g  (M  8°“ 

b  d 


(b  =  1,  3). 


Sy  (t)  =  g  (M  ^  ^  (b  =  1.  3). 

b  a  d 

Hg  Vgi7Q 

®(t)  =  g  (M  ^  ^)(s^(t),  (t))  Tg^ct)  5°'^. 

a  d 


Sy  (t)  =  g  (M  ^  )(s^(t),  0^  (t))  Tg^Ct) 

>'  d 

'^b'  ^B  *^B 

Sw  (t)  =  g  (M  ^  0G  ® 

b  a  °d 


Qg  '  Vg  q^ 

S„  "  (1)  =  g  (H  ^  S”'" 

b  ad 


(b  =  1,  ) 

a 


These  functions  are  the  contributions  from  the  body  gravity  force  to 
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the  total  quasi-coordinate  force  and  will  be  used  in  the  equations  of  motion 

"b 

(equation  C-16),  For  example,  the  simple  graph  S  will  be  included  on 

b 

the  right  side  of  the  V  equation  of  motion. 

O 


Suppose  that  F  is  the  vector  function  such  that  if  t  >  0,  F  (t)  is 
Hr  HP 

the  force  at  the  time  t  acting  on  the  point  HP  from  the  arresting  hook. 

The  point  HP  is  defined  in  Appendix  F.  Also,  at  the  time  t,  F  (t)  is  the 

H 

force  on  the  lateral  pivot  point,  LP,  from  the  lower  arresting  hook  segment 

and  is  directed  along  a  line  from  LP  to  HP.  F  (t)  is  equal  to  -  F  (t). 

H  HP 

Further,  suppose  that  r^p(t)  is  the  vector  function  such  that  if  t  0, 
rj^p(t)  is  the  vector  from  the  ground  reference  point  Q  to  the  point  labeled 
LP  at  the  time  t. 


The  virtual  work  done  by  the  arresting  force  is  therefore 


8 


W„(t) 


F„<t) 


(t). 


Now  suppose  that 


1  is  the  vector  function  such  that  if  t  >  0,  l„(t)  is  the  vector  from  the 
H  —  H 

jig  condition  location  of  the  airframe  reference  point  N  to  the  jig 
condition  location  of  the  point  labeled  H,  which  is  the  point  of  attachment 
of  the  arresting  hook  to  the  airframe,  at  the  time  t. 


1 


LP 


jig 


is  the  vector  function  such  that  if  t  >  0,  1  (t)  is 

L,r 

condition  location  of  the  point  H  to  the  point  LP  at 


the  vector 
the  time  t 


from  the 
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(t)o  is  the  vector  function  such  that  if  t  >  0,ct)„  (t)  is  the  bth 
airframe  modal  vector  at  the  time  t  for  the  point  labeled  H, 


"Vjf-  is  the  vector  function  such  that  if  t  ^  0,  (t)  is  the  bth 


airframe  modal  slope  vector  at  the  time  t  for  the  point  labeled  H  and 


I„  (t),  I„  (t)  and  I„  (t)  are  the  unit  vector  functions  such  that  if  t  >  0 
I„  (t)  has  the  direction  of  -  I, „(t),  I„  (t)  is  (t),  and 


I„  (t)  is  orthogonal  to  I„  (t)  and  (t)  and  is  directed  dovm 

H3 

with  respect  to  the  pilot.  See  Appendix  F  for  derivation  of  these  vectors 


Therefore,  the  vector  may  be  written  as 


Fj^pCt)  =  r(t)  +  Ij^(t)  + 

+  (t)  *  Ig  (t))  Ig  (t)  X  ^*8^^^^* 

Hb  3  3 


Based  on  the  assumption  of  small  airframe  deformations  the  vector 
function  r^p  may  be  differentiated  to  obtain 

Fj^p'Ct)  =  ?'(t)  +ng(t)  X  Ig(t)  +  (f2g(t)  •  Ig  (t))  Ig  (t)  X  I^p(t)  (11) 
(t)  q0‘’'(t)  +  (^g  (t)  •  Ig  (t))  Ig  (t)  X  I^^p(t)  qB‘''(t). 


The  vector  j’Lp'(t)  ®ay  be  used  as  a  replacement  for  Srg(t)  in  the 
virtual  work  equation  for  the  arresting  gear  force.  The  following 
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definitions  may  be  used  to  expand  the  vectors  in  equation  (11): 


F  ^  is  the  simple  graph  such  that  if  t  >  0, 
H 


F„(t)  =  F„  (t)  i,  at  the  time  t. 
n  n  D 


Each  of  and  c  =  1,  3;  d  =  1,  3  is  the  simple  graph  derived  in 

3  G 

Appendix  F  such  that  if  t  0, 

(t)  =  CLu°(t)  Tu^(t)  i  at  the  time  t. 

H3  d 


a^  and  are  the  simple  graphs  derived  in  Appendix  F  such  that  if  t  0, 
I  (t)  =  a  (t)  I_  (t)  +  b  (t)  I„  (t)  at  the  time  t. 


?'(t)  =  Vg°(t)  Ig  (t). 

c 

G 


Y„  (t)  =  Y  °  I  (t) 

"b 


It  follows  then  that  after  the  indicated  vector  operations  on  equation 
(11)  are  performed  the  following  coefficients  of  the  coordinate  velocities 
may  be  interpreted  as  the  quasi-coordinate  and  generalized  coordinate  force 
contributions  from  the  arresting  force: 
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B 


(b  =  1,  3) 


='‘>  =  TB^‘)  %df  S'"  S, 


(lB®(t)a„^t)  y  ^Ct)  Spg)  p/( 


B  '  '  H_ 
c  3 


xa„^(t)  7„^  (t)  (a^(t)  'yg*'(t)  +  b^(t)  yg''(t.))  § 


(c  =  1,  3) 


D 


^  ‘V  Tb®<«  a„V)  T„^t)  Sp,)  F„’(t) 

Hfa  c  3  h 


(t)  +  b^(t)  Tb''^')  S; 


(b  .  1,  N3g) 


Now  suppose  that  is  the  vector  function  such  that  if  t  0, 

F_  (t)  is  the  force  on  the  a  gear  wheel  from  the  ground  a.  the  time  t. 


Also,  suppose  that  r„  is  the  vector  function  such  that  if  t  >  0, 

^F 

a 

r„  (t)  is  the  vector  from  the  ground  reference  point  Q  to  the  point 

r 

^a  ® 

(the  a  wheel  tire  footprint  center  of  pressure). 


Therefore,  the  virtual  work  of  the  force  from  the  ground  on  the  tires 


C  q, 

'^W  (t)  =  X  Fp  ■  Or.  (t). 


a=1  a 
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It  is  noted  that  this  expression  is  the  virtual  woi k  done  by  the  ground 
forces  parallel  to  tiie  ground  plane.  These  ground  forces  are  derived  in 
Appendix  E.  The  generali::ed  forces  from  the  ground  which  are  normal  to  the 
ground  plane  will  be  derived  from  the  potential  energy  stored  in  the  tires. 


To  make  the  virtual  work  expression  above  more  usable  the  vector 

r„  (t)  must  be  expressed  in  terms  of  the  system  coordinates.  This  may 
a 

be  accomplished  as  follows: 


Suppose  w^  is  the  vector  function  such  that  if  t  0, 

w  (t)  is  the  vector  at  the  time  t  fr^m  the  point  A  (the  a  gear  axle 
U  2l 

a 

reference  point)  to  the  point  C  (a  point  that  is  common  to  the  a  gear  wheel 

a 

center'^lane  and  the  vertical  ground  force  vector) . 


Further,  suppos-'  that  v^  is  the  vector  function  such  that  if  t  0, 

V  (t)  is  the  vector  from  the  point  C  to  the  point  at  the  time  t.  In 
w  a 

a  j 

addition,  suppose  that  v^  is  the  simple  graph  such  tnat  if  t  0  then 

a 

Vy  (O  =  Vy  (t)  i  . 
a  a 


The  vector  r  (t)  ma/  be  expressed  as  follows; 
^F 

a 


r^  (t)  =  r(t)  +  (t)  +<$>gQ  (Sg(ty,  0^^  (t),  t)  (12) 

F  ^a  ^A  ■'a 


+  (t)  +  (s^(t),  t)  +  Cjlf,  (s^(t),  (t),  t)  ^(t)  +  w^  (t) 

a 


°A  ^ 
^b 


a 


If  the  vector  function  r  is  differentiated  and  it  is  assumed  that 

^F 
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the  magnitude  of  the  vector  w  (t)  is  independent  of  time  then 


?  ’(t)  =  Vg(t)  +  HgCt)  X  (t) 

a 

a 

“  a  a  a  a 

+  n,  (t)  X  (g^  (t)  +  (3g(t),  t)  +  (t)  +  (t)) 

3  a  a  a  a 

♦■  ^2  y  (  t  )  X  (  (  t  )  +  Vy  ( t ) ) 

a  a  a 

*'*80  <^3  ‘’'''  *’  ■’b'’’'*'*  *  ;3  =3'“’ 


(t),  t)  Qq  “’(t)  +  i  Vy  ’(t) 
a  a  a 


Since  at  the  time  t  the  vector  (t)  is  in  the  wheel  midplane  from 

a 

the  point  A  to  the  point  C  ,  it  follows  that  w-  may  be  expressed  by  a  unit 

a 

vector  times  the  ir^gnitude  of  w^  as  shown  by 

a 

(kG  X  I  )  X 

^2  ^2 

„  =  - : - : — ,  Wq  • 

a  |(k^  X  I^)  X  k^  I  a 

^2  ^2 


As  previously  defined,  the  unit  vector  k_  (t)  is  parallel  to  the  axle 

G 

®2 

of  the  a  gear  wheel  and  is  perpendicular  to  the  ground. 


Therefore  from  the  equation  for  k^  given  in  Appendix  G  it  follows 

®2 
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'*3  -  '’<^0  °  ‘  Tb"  I, 

^C  = -  ”c  ’ 

a  ^  a 


where 


^2  %  c  a^  a^  e 


O  b  ^  ^  7  3  /Q  f  e  y  2  2 

-*-  (Ag  '-^g  <5b  .-jq  ^g  Ob  ’ 

^2  S  c  a^  a^  e 


tQ  b  n  c  Y  3^  2,2  ,0.5 

+  (1  -  (^-Q  Qq  )  )  ]  . 

^2  ^b 


and  if  each  of  w^  ;  b  =  1 ,  3  is  a  simple  graph  such  that  if  t  >^  0,  then 
a 

«C  =  ”C  h 

a  a 


at  the  time  t. 


Thus,  if  r  '  is  used  for 
a 

ground  force  and  if  each  of  F„ 

U 


S  Tp  in  the  virtual  work  equation  for  the 

a 

is  a  simple  graph  such  that  if  t  >  0, 


a 

is  the  bth  component  of  the  a  gear  ground  force  such  that  at  the 


a 

time  t 


Fq  (t)  =  Fq  ‘’(t)  1^, 

a  a 

then  after  the  indicated  vector  operations  in  equation  (12)  are  performed, 
the  quasi-coordinate  and  generalized  coordinate  force  terms  for  the  ground 
force  are 
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B  G 

r<«  =  2.  Sbd  "g  “<'> 


(f  =  1,  3), 


a=1  f 


Sg  .('•>  =  2.  («.ah  h  "  S""  F.  ^t)  5 


c""  '  cdh 

a=  1  I  a 


Ga  -bg 


*  'ndh  Qg  <*A  *  *  S"'  F,  Nt)  S 


'cdh  ^G  '“A  ■'A  'a' 

fa  a  a 

e 


Ga  -bg 


*  ®fdb  "^B  ”c  ^G  ®f3b  "^B  ^G 


=G  -Qo  ^3h=  <«A  "  *  ^A  S””  F,  ^t)  S^g 


*  ^bOc  Tb"'''  “c'<«  "g 

a^  c  a  a 

*  ®d3b  ^G  ^W  ^G 

a^  c  a  a 

rz„ 

"g  ^G  Tb"^")  -C  "g 

a.  a  e  a  a 

2  p 

"  'o3b  Z^G  ’’'‘>  ^0  7b°'‘'  'w  <'^'  ■'g 

a„  a  e  a  a 

2  p 


B  G 

^G  ^G  "^B  ^^^^BG.  ^G  ^G  ^1 


a=1  a  e 

G 


(f  =  1.  Ng^), 


“a 

Sq  (t)  =  Qg  ^t)  Y/m  S,  '■•(ajt))  F^  ^t)  S^g. 

a  e  a  a  ® 

c 
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Sq  ^  f.(t)  =aQ  ®(t)  °(s^(t),  (9q  (t))  ‘'(t)  8 


bg 


(f  =  1,  Nqe  ), 
a 


where  b  ranges  from  1  to  2. 


To  derive  the  ground  force  normal  to  the  ground  plane  the  following 
functions  are  needed: 


3 

Suppose  that  is  the  simple  graph  such  that  if  t  >  0, 

G  — 

a 

F„  (x)  is  an  ordinate  of  F„  only  if  x  is  an  ordinate  of  the  simple  graph 

Cj  U 

a  a 

V,,  at  the  time  t  and  F.  ^(x)  is  the  normal  force  from  the  ground 
W  G  3 

a  a 

on  the  a  gear  tire  if  (t)  is  equal  to  x. 


Now  suppose  that  v  is  a  number  such  that  if  x  is  a  number  greater 

w 

*0 

than  or  equal  to  v„  then  F„  ^(x)  =  0.  Therefore,  if  u  is  a  simple  graph 

w  w  a 

^0 

such  that  the  point  (x,  u  (x))  belongs  to  u  only  if  x  is  an  ordinate  of  v 

a  a  w 

c 

and  X  is  less  than  and  u^(x)  is  the  potential  energy  stored  in  the  a 

^0 

gear  tire  corresponding  to  x  then 


\ 

/ 


=  -  •'x  '■g  ^  “I’ 

a 


or  alternatively 
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u^CVw  ](t)  =  -  J  F^,  dl. 


"'a 


and  since  the  generalized  force  associated  with  the  coo>'dinate  is 

a 

defined  by 


a 


it  follows  that 


V  ' 

^  ^<''W  =  -  ''g 

a  a  a 


The  other  two  components  of  the  ground  force  vector  are  derived  in 


Appendix  E. 


Suppose  that  is  the  velocity  of  the  wind  relative  to  the  ground. 
Further  suppose  that  Q  is  a  simple  graph  such  that  at  the  time  t,  Q(t)  is 
the  angle  of  attack  calculated  from 


Q(t)  =  tan 


^  (?'(t)  -  •  Ig^(t) 

(?'(t)  -  Vy)  • 


Also,  suppose  that  is  a  simple  graph  such  that  at  the  time  t,  jSit) 


is  the  side  slip  angle  calculated  from 


y8(t)  = 


(r'(t)  -  V  )  •  i  (t) 

-1  “  ®2 
j  - 

(?'(t)  -  •  I  (t) 
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Further  suppose  that  each  of  O  »  0  and  O  is  the  simple  graph  such 

that  ift>0,  j(t)is  the  elevator  rotation,  o  (t)  is  the  aileron 
—  e  a 

rotation  and  S  (t)  is  the  rudder  rotation  at  the  time  t. 
r 


-  12  3 
Now  if  F.  is  the  vector  function  such  that  if  (x„  ,  x„  ,  X-  ) 

A  D  O  D 

-  1  2  3 

is  in  V_  and  t  >  0,  ,  x^,  ,  x„  ,  t)  is  the  aerodynamic  force  per  unit 

12  3 

area  for  the  point  (x„  ,  x_  ,  x„  )  on  the  airframe  surface  at  the  time  t, 

D  O  D 

then  the  virtual  work  done  by  the  aerodynamic  force  over  the  surface  of  the 
airframe  S  is 

O 


s 


W^(t)  = 


F.(I  1,  12,  I  3,  t) 
B  B  B  B 


Sr 


R-(I  1,  12,  I  3,  t)  dS 
®  *B  *B  *B  ® 


Therefore,  if  q  is  a  simple  graph  such  that  the  time  t,  q(t)  is  the 
dynamic  pressure  and  if  P  is  the  air  density  then 

q(t)  =  1/2  P  Ir’(t)  - 

and  if  S  is  the  reference  area  for  the  aerodynamic  forces  and  if  it  is 
assumed  that  the  aerodynamic  force  is  independent  of  the  elastic  deformation 
of  the  airframe  then  there  is  a  vector  function  A  such  that 

=  q  s  (Xg\  Xg^,  Xg^,  a(t),/3(t),  8g(t),  SgCt),  S^(t), 

Ug\t),  N  tg^(t),  Ttg^Ct),  t). 


For  the  rigid  body  equations  of  motion  the  aerodynamic  forces  and 
moments  are  referenced  to  the  point  labeled  P  in  the  airframe.  Therefore, 
if  Ij^p  is  the  vector  function  such  that  if  t  0,  Ifjp(t)  is  the  vector  from 
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Che  jig  condition  location  of  the  reference  point  N  to  the 
aerodynamic  reference  point  P  at  the  time  t  and  if  is  the  vector 

rO 

1  2  3 

function  such  that  if  (x„  ,  x„  ,  x_  )  is  in  V  and  t  >  0, 

D  o  B  B  ““ 

12  3 

LnoCxn  ,  x„  ,  x„  ,  t)  is  the  vector  from  the  jig  condition  location  of  the 
rB  B  B  B 

12  3 

point  P  to  the  jig  condition  location  of  the  point  labeled  (x_  ,  x_  ,  x_  ) 

B  B  B 

at  the  time  t,  it  follows  that 

LgCXfi^  t)  =  r^p(t)  +  LpB(Xg^  Xg2,  t). 


This  expression  may  be  used  to  modify  equation  (5)  which  is  then  used 
to  replace  Or^  in  the  virtual  work  equation.  VHien  the  vector  operations 
and  integration  are  performed  as  indicated  in  the  virtual  work  equation  it 
is  seen  that  the  coefficients  of  the  velocity  terms  may  be  expressed  as 


b">  S,,  *  F/(t) 


“•b' 

^A  b('>  ■  Sbc 


(b  -  1,  3), 


(b  -  1,  3), 


(b  .  1,  Nj^). 


It  is  evident  then  that  the  aerodynamic  force  at  the  time  t  is 
F  ^(t)  i  (t)  and  the  aerodynamic  moment  at  the  time  t  is  I  (t)  H  '^(t). 

^  h  ^  \  ^ 

Further,  it  is  seen  that  (t)  is  the  generalized  force  at  the  time  t 
corresponding  to  the  bth  airframe  vibration  mode. 

Now  suppose  that  each  of  b=  1,  3  is  a  simple  surface  (aerodynamic 

force  coefficient)  such  that  the  aerodynamic  force  component  projection  at 
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the  time  t  on  i  (t)  may  be  expressed  by 
®b 

Fp'^Ct)  = 

q  S  Cp*'(a(t),  ^(t),  SgCt),  SgCt).  S^(t).  Qg^t),  Hg^Ct). 


(b  =  1,  3). 


Also,  suppose  that  each  of  C,,^;  b  =1,  3  is  a  simple  surface  (an 

aerodynamic  moment  coefficient)  and  each  of  b  =  1,  3  is  an  aerodynamic 

reference  length  such  that  the  aerodynamic  moment  projection  at  the  time  t 

on  i  (t)  may  be  expressed  by 
®b 

Mp^(t)  = 

q  S  I^'’  C^^(a(c),  S^(t),  S^ct),  S^(t),  Qg^t),  Qg^t),  .Qg\t)) 

(b  =  1,  3). 


Further  suppose  that  each  of  C*  ;  b  =  1,  N  is  a  simple  surface  (an 

Q 

aerodynamic  generalized  force  coefficient)  such  that 


Q^Nt)  - 

q  S  C^''(a(t),  /S(t),  Sg(t),  S3(t),  S^(t),  ng^(t),  ng^(t),  D  g^t)) 

(b  =  1.  Ngg). 


The  virtual  work  on  the  airframe  from  the  thrust  of  the  engine  is 
assumed  to  be  independent  of  the  airframe  deformation.  That  is,  it  is 
assumed  that  if  TL  is  a  point  on  the  thrust  line  of  the  engine(s)  then  the 
airframe  deformation  in  the  direction  of  the  thrust  force  is  zero. 
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Suppose  that  is  the  vector  function  such  that  if  t  ^  0, 

F^(t)  is  the  thrust  force  vector  at  the  time  t  and  the  magnitude  of  F^(t) 
is  a  number.  The  thrust  force  vector  P.j(t)  may  be  expressed  in  component 
form  as 

F.j.(t)  = 

with  F„(t)  •  (t)  =  0. 

T 

Also  suppose  that  is  the  vector  function  such  that  if  t  0,  r.j,(t) 
is  the  vector  from  the  Jig  condition  location  of  the  airframe  reference 
point  N  to  the  jig  condition  location  of  the  point  TL  a  the  time  t. 
Therefore,  if  r^^^  is  the  vector  function  such  that  if  t  0, 
vector  from  the  ground  reference  point  Q  to  the  point  TL  at  the  time  t,  then 
the  virtual  work  done  by  the  thrust  force  is 

Sw^(t)  =  F.j(t)  •  5?^^(t). 

The  vector  written  in  expanded  form  excluding  the 

deformation  terms  as 

rTL(t)  =  r(t)  +  I^(t). 

Thus , 

rTL’(t)  =  r'(t)  +^^g(t)  x  i^(t). 

The  vector  r  '(t)  may  be  used  as  a  replacement  for  Sr_,  (t)  in  the 
X  L«  TL 

virtual  work  equation.  It  follows  then  that  if  in  component  form 

iT(t)  =  I  (t). 
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the  quasi-coordinate  force  contributions  from  the  thrust  force  are 


S 


be 


o  —  p  ^  p  1  ^ 

T  2  "  T  2bc  T  • 


(b  =  1,  3), 


Now  suppose  that  i^  is  the  vector  function  such  that  if  t  0, 

^3 

ig  (t)  is  the  unit  vector  at  the  time  t  which  is  parallel  to  the  a  gear 

shock  strut  axis  of  symmetry  and  directed  downward  with  respect  to  the 

piston  as  shown  in  Figure  B-2  of  Appendix  B.  Further,  suppose  that 

F  is  the  vector  function  such  that  if  t  >  0,  F-  (t)  is  the  shock  strut 
a  a 

force  (see  Appendix  B  for  derivation  of  the  shock  strut  force)  at  the  time 

t.  Also,  suppose  that  Fg  is  the  simple  graph  such  that 


Therefore,  the  selection  of  the  a  gear  stroke,  s  ,  as  a  generalized 
coordinate  permits  the  virtual  work  from  the  shock  strut  force  to  be 
expressed  by 

Swg(t)  = 


"G 

1 


( t)  Os  ( t) . 


Therefore,  the  contribution  from  the  shock  strut  force  to  the 


generalized  force  in  the  a  gear  equation  of  motion  is 

3  * 
a 

Sg  (t)  =  Fg  (t). 

a 


44 


The  deformation  of  the  airframe  and  the  deformation  of  the  landing 
gears  are  assumed  to  be  represented  by  normal  vibration  modes. 

Consequently,  the  stiffness  matrices  for  the  airframe  and  landing  gear 
deformation  equations  of  motion  are  diagonal.  For  simplicity,  it  is  assumed 
that  the  damping  is  viscous  and  the  damping  matrices  for  these  equations  are 
diagonal.  Based  on  these  assumptions  the  following  definitions  are  made. 

Suppose  that  b  is  an  integer  in  [1,  and  that  the  number  is  the 

generalized  stiffness  corresponding  to  the  bth  airframe  modal  equation. 
Further  suppose  that  the  number  C„^  is  the  generalized  damping  corresponding 

O 

to  the  bth  airframe  modal  equation. 

Therefore,  the  generalized  force  from  the  airframe  modal  stiffness  and 
damping  for  the  bth  mode  is 


Also,  suppose  that  c  is  a  positive  integer  in  [1,  N  ] ,  is  the 

Ob  C 

a  a 

fully  compressed  stroke  of  the  a  gear  and  K  is  a  simple  surface  such 

G 

a 

c 

that  if  X  is  in  [0,  s„  ] ,  y  is  in  [-TT,  TT  ) ,  K  (x,  y)  is  the  generalized 

G  G 

a  a 

c 

stiffness  for  the  cth  mode  of  the  a  gear  corresponding  to  the  a  gear  stroke 

X  and  castor  angle  y.  Further,  suppose  that  C  is  a  simple  surface  such 

G 

a 

c 

that  if  X  is  in  [0,  s^  ] ,  y  is  in  [-TT,  TT],  (x,  y)  is  the  generalized 

a  a 

c 

damping  for  the  cth  mode  of  the  a  gear  corresponding  to  Che  a  gear  stroke  x 
and  castor  angle  y. 
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The  generalized  force  from  the  a  gear  modal  stiffness  and  damping  for 
the  cth  mode  is 


GF  c 


(t)  =  -  K-  (s  (t).  Or  (t))  dr  ‘^(t)  S, 


G  '  a 
a 

G 


-  '^0  <V‘>-  ^0  Sed- 

a  a  a 

c 


As  stated  previously,  it  is  usual  for  only  the  nose  gear  to  be  capable 
of  castoring.  However,  the  final  equations  may  be  written  more  concisely  if 
all  the  gears  are  assumed  to  have  this  degree  of  freedom. 


The  gear  castoring  moment  is,  in  general,  a  nonlinear  function  of  the 

castor  angle  function  for  the  a  gear,  0  ,  and  its  derivative.  The 

G 

a 

specific  functional  form  of  this  moment  is  usually  configuration  dependent, 

Therefore,  to  simplify  the  current  derivation  it  is  assumed  that  is  a 

simple  surface  such  that  if  t  >  0,  M_(^_  it) ,  O,  „  (t))  is  the  generalized 

—  u  Cj  G 

a  a 

moment  associated  with  the  generalized  coordinate  Or-  time  t. 

G 

a 


Thus,  the  castoring  generalized  force  may  be  written  as 


^(t)  =  -  (t),  Hp.  (t)). 
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The  Constraint  Condition 


In  the  previous  section  the  coordinate  v,,  was  introduced  when  the 

W 

a 

ground  forces  acting  on  the  system  were  defined.  As  a  consequence  there  is 
a  redundancy  in  the  coordinates  specifying  the  position  of  the  system.  This 
redundancy  may  be  easily  handled  by  use  of  a  constraint  condition  and  a 
Lagrangian  multiplier  as  illustrated  in  Appendix  C. 


The  constraint  condition  may  be  stated  as  follows:  If  at  time  t  the  a 
gear  is  on  the  ground  then  the  component  of  the  velocity  r  '(t),  from 


equation  (13),  normal  to  the  ground  is  zero.  Therefore, 


r-  ’(t)  ’  I,  =  0. 
a 

As  seen  from  equation  (13)  this  constraint  condition  has  the  form  of 
equation  (C-13)  in  Appendix  C.  Consequently,  as  seen  from  equation  (C-14) 
the  constraint  function  may  be  written  in  the  following  form: 


’^B  ’^B^’^'G 


BE. 


=  0. 


After  the  indicated  vector  operations  have  been  performed,  the  terms  to 
be  used  on  the  right  side  of  equation  (C-16)  are 
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Equations  of  Motion 


All  of  the  input  terms  -o  the  equations  of  motion  (equation  (C-l6)  of 
Appendix  C)  have  now  been  derived.  The  definitions 


"b  "b  '^b  "b  "b  ''b 

^  b  "  b  b  b  b  b 

Hb  ^B  ^B  ^B 

^  b  =  b  b  ■"  b  ■"  ^A  b  b 

r^TQ  Cl  Q  Cl  Q 

^  \  ■*■  ^G  ■"  ^GC 


Cl y  Cl  y 

s  ^  -  s  ^ 


s'  s  *  s'  s 
a  a  a  a 

S  =  S.,  +  +  S 


W 


S  ’ 


‘^b'  '^b'  ‘^B*  ‘^b 


'  ^b’  '^b’ 


^  b  =  b  ■"  b  ^G  b  ■"  ^A  b  ■"  ^BF  b 


‘^G  ’  '^G  ’  ^G  ’  '^G  ’ 

s®  -S  ^  +S  ^  +'*  ^ 

^  b  ■  b  ^  ^G  b  *  “GF  b 


(b  =  1,  3). 

(b  =  1,  3), 


(b  ,  1,  Njg), 

(b  ■  1,  ), 


permit  the  equations  of  motion  to  be  expressed  in  the  following  form: 


0  =  S, 


.X 


a  a;v. 


(K.y  1)'  K.„  2  K.y  3  = 


'P  ;V3-  '•'B  -.Vg 


1 
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2 


V  N 
''b  G 


(K.y  2)'  +^2^  ^  ^  ^ 


V  N 
B  G 


(K.  3)'  -O.^  It.V  '  "iVg^  '  ^  3 


' ''»'  -^b' 


"g 


B  1  *  ?,  ^a  ^a-.fto' 


'Sa„^''  '  ''b'  -  ''b’  Sv/  ^Ob'  -^b’  = 


^G 


*  'b'  Sv  b  -fig^  Sn ;  -  Ob’  = 


'^G  ,  . 

s  3  +  X  A  F.q  3. 

^  a=1  ’'^‘b 


'^n.  '■  -  ’'^1.  = "  '  *  ^n, 


'G  ‘•G 

a  a 


(a  =  1,  Nq) 


(K.o  )'  =  S 

>0  6y 

a 


(a  =  1,  Nq) 


(K,3^,.'  -S^  =  S 


(a  =  1,  N^) 
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(K  b,)'  -  b  =  S  K  F.^  b, 

.n  .  .n  b  a  ;qg 


‘B 


(b  =  1,  Ngg)^ 


(K.^  b,)'  -  K.^  b  =  S 


•^G 


;<3r 


■"  (a  =  1,  and  b  =  1,  Ngg  ). 


b  a  a;q 


These  equations  may  be  rewritten  by  use  of  the  expanded  expression  for 
the  kinetic  energy  and  the  generalized  and  quasi-coordinate  forces. 


The  first  equation  of  motion  above  permits  the  calculation  of  the 
Lagrangian  multiplier  X_.  Since 

cl 


a;v^ 


=  1.0 


and 


a 


‘^G 

then 


=  -  F 


G  ’ 
a 


a 

Thus,  the  Lagrangian  multiplier  X  is  seen  to  be  the  component  of  the 
force  on  the  a  gear  tire  from  the  ground  that  is  normal  to  the  ground. 
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The  V„  equation  is 

D 


V  V 

B 


^G  '^B  ^  B 


^G  '^B^B 


1  o\-s^  ^a'  ^b""  1  o'>;0^  ^G^  ^B 


^G  ^B  ^G 


*^G  '^B  ^  G 


Z  1  ^G  ^  ^G  '  •*■  2  I  ^;s 

azi  a  a  a=1  a 


o  1  ^  / 

^  •*  r* 

SI  G 


^G  ^B  ^  G 


(0.  )^ 


S,  '  >^00  ‘"“a 


"g  "b  =a 


"a  "b  "a 


*  I  (M  ,  )[3^,  (9<,  1  a,,"  •►  I  (M  ,  ).  (a.') 


"g  ''b  =a 


'b  ■’b 


X  <"  ,  >;f  «G  =a'  *  '«  1  c>'=.’  ^G  ' 

azi  ’'^G  a  a 


Nq]  Qb  " 


*^G  '^B  '^B 

I,  1  o);s  ^b"’ 

a=l  a 


«G  '’b 


S,  ‘”  '  <=^-do  ^°a  ' 


'^G  '^B  ^G 


*^G  ^B  '^G 


*  2  <«  ,  %)t=a.  &,  1  lo  °"  *  1  <M  ,  ‘c*; 

azi  a  a  a=l 


3  '  ' 

3  a  G 


"g  "b  '^G 


+  X  ^  qK  C  ^Q  ^G  * 
a=1  '  '  ^G  a  a 
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Vg  Vg  Vg  Hg 

((M  ,  ,)  Vb^  *  (M  J  ^  »o’^B 

a 

Nq  Vg  Hq 

*S<«3  ^[3^,  0b>^G 

a=i  a  a 


"o  ’'b  \  ^  Vb  Qb 

*  X  (M  2  )[3,.  6^  1  3^'  *  (M  B  ^)[3  (9 

*  a 


)C3^,  9b  i  3  =  1.  Nb)  Qb”' 
a 


»o  ''b  Og 

2  ^0  ’  Oo  °'> 

a=  1  a  a 


Vg  Vg  Vg  rig 

1  V  (M  3  ^)[s^,  ;  a  =  1,  N^]  Hg® 

a 

^  G 


S  (M  ,  ^[s  .  Or.  la 


a»  <^Q  j  ufcQ 

a  a 


Nq  Vg 

^  S  (M  )[s  0  ]  3  . 

a=1  ^  ^  a 


'^B  '^B 


3  c^'^®a’  6^g  5  a  =  ^  Nq]  qg°' 


Nq  Vg 

2  («  3  %)C3^.  9b  J 

a=i  a  a 


V  V 
B  B 


*  (  Tb^  Sbb)  h.  8,b. 
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The  Vg  equation  is 


■b  h 

(M  ,  ,)  +  (M  2  ;  a  =  1,  N^]  Hg®' 


1  r  B 


(15) 
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The  V_^  equation  is 
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Initial  Conditions 


The  aircraft  is  assumed  to  be  oriented  at  t  =  0  such  that  one  point  of 
the  aircraft  has  just  made  contact  with  the  ground  surface.  The  candidate 
points  for  contact  are  the  landing  gear  tires  and  the  arresting  hook.  There 
are  some  cases  when  there  are  multiple  points  of  simultaneous  contact  at  the 
time  t  =  0. 

The  transformation  ^^*^(0);  b  =  1,  3  and  c  =  1,  3  may  be  established 
from  the  Euler  angles  ^(0),  ^(0)  and  cj^iO)  which  are  the  initial  yaw, 
pitch  and  roll  angles  of  the  airframe.  Appendix  G  defines  the  elements  of 
this  transformation. 

The  rigid  body  translational  velocity  of  the  airframe  reference  point  N 
at  t  =  0  is 

?'(0)  =  V„°(0)  (0)  =  C^’(O)  I  =  d'(0)  i,  +  3'(0)  +  v'(0)  I,. 

o  D  b  1  2  3 

c 

Normally,  the  horizontal  velocity,  the  vertical  velocity  or  sinking 
velocity  and  the  lateral  velocity  of  the  point  N  with  respect  to  the  ground 
surface  (i.e.  ^^'(0)  ^^’(0)  and  ^^'(0)  are  specified  at  t  =  0, 

Therefore,  the  body  axis  components  of  the  initial  velocity  are 

d 

Also,  when  the  initial  angular  velocity  is  not  zero,  the  components  of 
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the  initial  angular  velocity  may  be  expressed  in  terms  of  the  initial  Euler 
angle  rates  as  follows: 

Hq^O)  =<^'(0)  -  sin(0(O))  ^'(0). 

=  cos(cf>(0))  0<(O)  +  003(0(0))  sin(C^(0))  V^«(0). 

(0)  =  -  sin(</:>(0))  0'(O)  +  cos(0(O))  cos((^(0))  ^'(0). 

The  following  bth  airframe  modal  displacement  for  b  in  [1,  Ngg]  at 
t  =  0  is  derived  from  equation  (20).  It  is  assumed  that  there  may  be  cases 
where  the  initial  castor  angle  is  not  zero.  It  is  further  assumed  that  the 
functions  that  are  dependent  on  the  stroke  of  the  landing  gears  can  oe 
evaluated  with  the  initial  stroke  equal  to  zero  then 

Vb  % 

(g  (M  ^  ^j)(0,  0Q  (0);  d  =  1,  Nq)  8°®  ♦  Q^^). 

B.  d  e 

b 

The  initial  position  of  the  point  N  with  respect  to  the  ground 
reference  point  Q  is  defined  as  follows: 

5'<0)  =  <!«»  =  ?(»)  •  I,- 

^^(0)  =  s(0)  .  r(0)  •  ij. 

^3(0)  =  v(0)  =  r(0)  ■ 

The  number  d(0)  is  normally  zero,  the  number  3(0)  is  used  to  orient  the 
point  N  laterally  with  respect  to  the  point  S  (see  Appendix  F)  and  the 
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nufflber  v(0)  may  be  derived  from  equation  (12)  for  the  case  where  the  a  gear 
tire  has  just  contacted  the  ground  surface.  It  is  found  from  equation  (12) 


that 


v(0) 


- 


-4^bo  “ 

‘  a  a 

c 


7  5(0  q„^0) 


G  7b5(0) 

^b  ° 


-^r  '"(0-  9r  (0))  CL.  "*(0)  T.^(O)  q  ‘'(O) 


"c 

a 


t  = 

v(0) 


The  equation  for  v(0)  may  be  obtained  a'  follows  for  the  case  whe’^e  at 
0  the  arresting  hook  i*''"?  just  n.dde  contact  with  the  ground  surface; 


7b^(0,  -  < 
b 


HV 


3in(^  )  y^ho) 


.  1„,  003 


>  Te 


’(0)) 


The  a  gear  stroke  at  t  =  0  is  determined  from  the  following  equation 

which  is  derived  from  equation  (21)  and  the  equation  for  the  number  F  (0) 

a 

from  i'.ppendix  B: 


^  b  c^s  Nq)  i^g'^O)  Hg^^LO) 

’  a  d 

^B^B 

"  b  c’|3  6’g  <“'■  “  '  '•  "g’ 


B  a  ^ 

*  8  (M  ^  )(0,  (0))  yg^CO) 
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The  hydraulic  fluid  expansion  in  the  a  gear  fluid  chamber  and  the  a 
gear  snubber  chamber  may  te  obtained  at  t  =  0  from  Appendix  B  as  follows 


/S.CO)  =  - 


p,  U  W)  v„  (0) 

a  a 


■Si 


(0)  =  - 


Pa  (3  (0))  V  (0) 
a  S 

a 


The  following  a  gear  bth  modal  displacement  for  b  in  [ 1 ,  ]  at  t 

GE 

a 

is  deri/ed  from  equation  (23); 


^0) 


(g  (h  \)(o,  0  (0))  y  gcd) 
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I-H  t-H  fO  I  1^  (M  In  I  H 


Summary  of  Vector  Function  Relations 


rj  8  d 


7S 


r 

I 


a 


w 

a 


b 


B  ’  ^  B  ^B 
b  bo 

BG  ’  ^  BG  ^G 
a.  a 
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SIMPLIFICATION  OF  EQUATIONS 


The  equations  of  motion  (equations  (14)  through  (24))  in  combination 
with  the  equations  developed  in  the  Appendices  provide  considerable 
flexibility  for  calculcting  the  loads  associated  with  landing  impact,  A 
minimum  of  assumptions  were  made  during  the  course  of  the  development  of 
these  equations  to  help  ensure  this  flexibility.  In  many  cases,  however, 
the  equations  may  be  simplified  by  virtue  of  the  gear  design.  For  example, 
for  the  case  where  the  nose  gear  power  steering  unit  is  engaged  and  is  able 
to  maintain  the  castor  angle  at  zero,  then  considerable  simplification  may 
be  realized. 

There  are  several  simplifications  that  may  he  made  to  the  equations  of 
motion  if  the  .-^ears  are  not  articulated  (i.e.  there  is  a  one  to  one 
correlation  of  the  shock  strut  stroke  and  axle  motion) .  These 
simplifications  are  readily  determined  from  an  .wamination  of  the  equations 
of  motion. 

There  are  other  simplifications  that  may  be  realized.  For  example,  it 
is  usually  permissible  to  assume  that  the  airframe  vibration  modal  effective 
masses  are  independent  of  the  shock  strut  stroke.  Therefore,  when  this 
approximation  is  used  in  conjunction  with  airframe  normal  modes  that  are 
orthogonal  to  the  rigid  body  motion,  the  complexity  of  the  equations  of 
moofon  is  reduced. 

Further,  it  may  be  possible  to  delete  other  terms  by  computing  the  mass 
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terms  in  Appendix  D  ana  estimating  the  contribution  of  the  terms  involving 
these  masses  with  approximations  of  the  multiplying  functions  which  are 
found  in  the  equations  of  motion. 


82 


REFERENCES 


1.  Keith,  J.S.,  "Methods  in  Structural  Dynamics  for  Thin  Shell  Clustered 
Launch  Vehicles,"  FDL-TDR-64-105,  April,  1965 

2.  Whittaker,  E.T.  "A  Treatise  on  the  Analytical  Dynamics  of  Particles  and 
Rigid  Bodies,"  Cambridge  University  Press,  1961 

3.  Wall,  H.S.,  "Creative  Mathematics,"  University  of  Texas  Press,  1963 

4.  Craig,  H.V.,  "Vector  and  Tensor  Analysis,"  McGraw  Hill  Book  Co.,  1943 

5.  Goldstein,  H.,  "Classical  Mechanics,"  Addison-Wesley  Publishing  Co.,  1959 


1 

i 


83 


APPENDIX  A 


MATHEMATICAL  NOTATION 


In  the  text  of  the  report  there  is  frequent  use  of  the  terms  simple 
graph  and  simple  surface.  This  is  done  to  emphasize  their  functional  and 
geometrical  significance. 

The  statement  that  "f  is  a  simple  graph"  means  that  f  is  a  point  set  no 
two  members  of  which  have  the  same  abscissa.  The  number  f(x)  is  the 
ordinate  of  that  point  of  the  simple  graph  f  whose  abscissa  is  x. 

The  statement  that  "T  is  a  transformation"  means  that  T  is  a  collection 
of  one  or  more  ordered  pairs  no  two  of  which  have  the  same  first  term. 

The  statement  that  "F  is  a  simple  surface"  means  that  F  is  a 
transformation  from  a  point  set  to  a  number  set.  If  the  point  P  =  (x,  y) 
and  the  number  z  constitute  a  point  of  F  then  z  =  F(x,  y) . 

The  use  of  the  bracket  product  notation  has  been  used  extensively  in 
this  report.  The  following  definition  of  a  bracket  product  may  be  used  to 
illustrate  the  essential  features  of  this  important  functional  description. 

The  statement  that  "the  simple  graph  F[g^ ,  g^J  is  the  bracket  product 
of  F  of  g^  and  g^  (F  is  a  simple  surface  and  each  of  g^  and  g^  is  a  simple 
graph)"  means  that  there  is  a  number  x  such  that  (g.|(x),  g^Cx))  is  in  the 
x,y  projection  of  F  and  If  x  is  such  a  number, 
fLg^.  g2Kx)  =  F(g^(x),  g^Cx)). 
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An  extremely  useful  concept  in  the  handling  of  functions  and  their 
derivatives  is  the  identity  graph.  This  graph  contains  the  point  (0,0)  and 
has  slope  equal  to  one.  Consequently  the  statement  that  I  is  the  identity 
graph  means  that  if  x  is  a  number  then  I(x)  =  x.  The  identity  graph  I 
appears  in  the  text  with  a  subscript  to  provide  clarity  when  used  in 
combination  with  simple  surfaces. 

With  the  help  of  the  identity  graph  the  definition  of  the  derivative  of 
a  simple  surface  may  be  stated  as  follows:  The  statement  that  is  the 
1-derivative  of  the  simple  surface  F  means  that  F^'  is  the  simple  surface  to 
which  the  point  ((x,  y) ,  2)  belongs  only  if 
z  =  {F[I,  y]}'(x). 

In  addition  the  statement  that  F^'  is  the  2-derivative  of  the  simple  surface 
F  means  that  F^’  is  the  simple  surface  to  which  ((x,  y) ,  2)  belongs  only  if 
z  =  {F[x ,  I] } ' (y) . 


These  definitions  lead  to  the  notation  used  in  the  text  of  the  report 
for  partial  derivatives.  To  illustrate  this  suppose  that  each  of  f,  g^,  and 
g^  is  a  simple  graph  and  F  is  a  simple  surface  such  that 
f  =  F[g^ .  ^2^' 

The  statement  that  F_  is  the  partial  derivative  of  F  with  respect  to  g 

>  8  1 

means  that 


Further,  the  statement  that  F_  is  the  partial  derivative  of  F  with  respect 


f  i^2 


to  g^  means  that 
^:g2  ^  ®2^- 
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The  identity  graph  is  also  useful  in  the  integration  process.  In  the 
general  case  f  is  a  simple  graph  bounded  on  the  interval  [a,  b]  and  the 
simple  graph  g  is  nondecreasing  on  r.a,b]  then  f  is  g  integrable  on  [a,b] 
(Reference  (3)).  This  is  denoted  by 


In  the  special  case  where  g  =  I  (the  identity  graph)  this  integral 
(i.e.  the  integral  of  f  with  respect  to  I)  is  the  well  known  Reimann 
integral . 

Vectors  are  denoted  by  a  bar  over  the  symbol.  When  the  vector  is 
expressible  in  terms  of  a  unit  vector  set,  for  example,  I  ;  a  =  1,  3  such 

3 

that  if  F  is  such  a  vector  and 
F  =  F  I^  (sum  on  a  is  implied), 

3  •• 

then  each  of  F  ;  a  =  1 ,  3  is  referred  to  as  a  component  of  the  vector  F. 


Frequent  use  is  made  of  the  "e  system”  for  representing  vector 
operations.  This  system  is  discussed  in  detail  in  Reference  (4).  The  e 
system  may  be  defined  in  the  following  manner: 

If  each  of  r,  s  and  t  is  an  integer  in  the  interval  [1,  3]  then  e  is  -1, 

1  s  t> 

0  or  +1  according  to  the  formula 

=  1/2  (r  -  s)  (s  -  t)  (t  -  r) , 

and 
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It  is  seen  that  is  skew-symmetric  in  each  pair  of  its  Indices. 


The  e  system  is  used  in  the  main  body  of  the  report  to  express  vector 
cross  products  in  terms  of  vector  components.  For  example,  if 

A  =  A®  I 

s 

and 

S  =  b‘ 

then 

A  X  B  .  i  =  e  ^  A®  B^. 
r  rst 


The  Kronecker  deltas  used  in  the  text  are  defined  as  follows: 


5 


abc 

rst 


'rst 


abc 

e 


Sab  _  ^  abc 
rs  “  ^rsc' 


Also, 

S*”  -  5,.  =  1 
8*^  =  8,,  =  0 


if  a  is  equal  to  b  and 


if  a  is  unequal  to  b. 
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APPENDIX  B 


DERIVATION  OF  THE  SHOCK  STRUT  FORCE  EQUATIONS 

The  following  shock  strut  characteristics  are  assumed  for  the 
derivation  of  the  equations  for  the  shock  strut  pressure  and  force  functions 
in  terms  of  the  quasi-coordinate  velocity  component  functions  of  the 
airframe  and  the  generalized  coordinate  displacement  functions  for  the 
landing  gear: 

(1)  The  shock  strut  contains  a  hydraulic  fluid  and  a  gas. 

(2)  During  shock  strut  compression  energy  is  stored  in  the  gas. 

(3)  During  shock  strut  compression  and  extension  energy  is  dissipated 
in  the  hydraulic  fluid  and  the  coulomb  friction  forces  from  the  shock  strut 
bearings . 

(4)  The  shock  strut  is  equipped  with  an  orifice  for  control  of  the 
hydraulic  fluid  flow  during  compression  and  the  area  of  the  orifice  is 
assumed  to  be  dependent  only  on  the  stroke  of  the  shock  strut. 

(5)  The  shock  strut  is  equipped  with  a  snubbing  chamber  (see  Figure 
B-1)  for  energy  dissipation  during  extension  of  the  shock  strut  and  the 
orifice  area  of  the  snubbing  chamber  is  controlled  by  a  pressure  relief 
valve . 

(6)  The  shock  strut  has  a  single  stage  air  chamber. 

(7)  The  hydraulic  fluid  is  compressible  and  the  bulk  modulus  of  the 
hydraulic  fluid  is  assumed  to  be  a  number. 

To  provide  some  additional  flexibility  for  modeling  potential  shock 
strut  configurations,  two  shock  strut  geometries  that  have  the 
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characteristics  described  above  will  be  included  in  the  derivation. 


The  terms  used  in  the  shock  strut  force  derivation  are  defined  below. 
Figures  B-1  and  B-3  may  be  referred  to  for  illustration  of  the  geometrical 
quantities.  The  required  numbers  are  defined  as  follows: 


is  the  ambient  pressure. 

AnB 


A  is  the  hydraulic  area  of  the  a  gear  fluid  chamber. 
H 

a 


is  the  net  area  exposed  to  the  gas  pressure  in  the  a  gear. 


Aq  is  the  area  of  the  hole  in  the  main  orifice  plate  in  the  a  gear. 


A^  is  the  hydraulic  area  of  the  gear  snubbing  chamber  in  the  a  gear. 
5 

a 

K  is  the  bulk  modulus  of  the  hydraulic  fluid, 
n 


is  the  mass  density  of  the  hydraulic  fluid. 


V  is  the  gas  volume  in  the  a  gear  gas  chamber  with  the  shock  strut  fully 

a 

extended . 


V  is  the  volume  in  the  a  gear  fluid  chamber  with  the  with  the  shock  strut 
a 

fully  extended. 
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K_„  is  the  spring  rate  of  the  a  gear  piston-cylinder  interface  when  the 
Bu 

a 

strut  is  fully  extended. 


is  the  damping  rate  (assumed  viscous)  of  the  a  gear  piston-cylinder 

bu 

a 

interface  when  the  strut  is  fully  extended. 


V„-  is  the  fluid  volume  in  the  a  gear  shock  strut, 
r  L 

a 


is  the  main  orifice  coefficient  for  the  a  gear, 
a 


Cp  is  the  snubbing  chamber  orifice  coefficient  for  the  a  gear. 

s 

a 

is  the  mass  of  the  rotating  parts  on  the  a  gear  axle, 
a 


Sj,  is  the  stroke  of  the  a  gear  when  fully  compressed. 


is  the  upper  shock  strut  bearing  friction  coefficient. 


is  the  lower  shock  strut  bearing  friction  coefficient. 


The  vector  functions  needed  to  define  the  shock  strut  force  are  defined 
as  follows; 


ig  is  the  unit  vector  function  such  that  if  t  >_  0,  i^  (t)  is  the 

unit  vector  which  at  the  time  t  is  parallel  to  the  shock  strut  axis  of 
symmetry  and  directed  downward  with  respect  to  the  piston  as  shown  in  Figure 
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B-2. 


j 

I 

F„  is  the  vector  function  such  that  if  t  >  0,  F  (t)  is  the  ohock 

a  a 

strut  force  on  the  a  gear  piston  at  the  time  t. 


F  is  the  simple  graph  such  that  if  t  0, 
a 

Fg  (t)  =  Fg  (t)  (t)  at  the  time  t. 

a  a  a^ 

v„  is  the  vector  function  such  that  if  t  >  0,  v„  (t)  is  the  velocity  of  the 

H  —  M 

a  a 

hydraulic  fluid  passing  through  the  a  gear  orifice  at  the  time  t. 


v„  is  the  simple  graph  such  that  if  t  >  0, 
n  — 

a 


v„  (t)  =  -  Vu  (t)  i-  (t)  at  the  time  t. 
n  no 

a  a  a^ 


V,,  is  the  vector  function  such  that  if  t  >  0,  v„  (t) 

a  a 

is  the  velocity  of  the  hydraulic  fluid  passing  through  the  a  gear  snubbing 

orifice  at  the  time  t. 


V  is  the  simple  graph  such  that  if  t  >  0, 

“s 

a 

V  (t)  =  -  V  (t)  i  (t)  at  the  time  t. 

“s  a, 

a  a  3 


F„  13  the  vector  function  such  that  if  t  >  0,  F„  (t)  is  the  shock 
F  —  F 

a  a 

strut  friction  force  on  the  a  gear  piston  at  the  time  t. 


P  is  the  simple  graph  sucn  that  if  t  >  0 


Fp  (t)  =  -  Fp  (t)  (t). 

a  a  a,^ 
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F  is  the  vector  function  such  that  if  t  >  0,  (t)  is  the  a  gear 

a  a 

axle  force  from  the  axle-wheel  interface  at  the  time  t. 


Each  of  F^  :  b  =  1 ,  3  is  a  simple  graph  such  that  if  t  >  0, 

F„  (t)  =  F-  (t)  (t)  at  the  time  t. 

G.  G,  vj 

A  A  a 

a  a.  c 

b 

M-  is  the  vector  function  such  that  if  t  >  0,  (t)  is  the  a  gear 

K  ~  \ 

a  a 

axle  moment  from  the  axle-wheel  interface  at  the  time  t. 


and  each  of  M-  ;  b  =  1 ,  3  is  a  simple  graph  such  that  if  t  >  0, 

A  \ 

A,  _ 

Mq  (t)  =  Mq  (t)  iy  (t) 
a  a^  1 

=  Mq  (t)  3q  (t) 

A  a 

The  additional  simple  graphs  required  to  determine  the  shock  strut 
force  are  defined  below: 


3^  is  the  simple  graph  such  that  if  t  0,  s^{t)  is  the  a  gear  stroke  at  the 
time  t. 


Ap  is  the  simple  graph  such  that  if  x  is  in  [0,  s^  ],  Ap  (x)  is  the 
a  a 

area  of  the  metering  pin  at  the  a  gear  stroke  x.  Ap  (x)  is  the  cross 


sectional  area  in  the  plane  where  the  orifice  plate  hole  area  A^  is 

a 
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measured . 


A.,  is  the  simple  graph  such  that  if  x  is  in  [0,  s  ],  A  (x)  is 
a  a  a 

the  a  gear  net  orifice  area  of  the  main  orifice  at  the  a  gear  stroke  x. 


A.,  is  the  simple  graph  such  that  if  x  is  a  number,  A„  (x)  is  the 

h  S 

a  a 

a  gear  net  orifice  area  for  anubbiug  at  the  snubber  orifice  differential 
pressure  x. 


is  the  simple  graph  such  that  if  t  0,  volumetric 

expansion  of  the  hydraulic  ^luid  in  the  a  gear  fluid  chamber  at  the  time  t. 


is  the  simple  graph  such  that  if  t  0,  (t)  is  the  volumetric  expansion 

^a  a 

of  the  hydraulic  fluid  in  the  a  gear  snubbing  chamber  at  the  time  t. 


V„  is  the  simple  graph  such  that  if  t  >  0,  V  (t)  is  the  a  gear  fluid 
a  a 

chamber  volume  at  the  time  t. 


V  is  the  simple  graph  such  that  if  t  >  0,  V„  (t)  is  the  a  gear  snubbing 

a  a 

chamber  volume  at  the  time  t. 


p  is  the  simple  graph  such  that  if  x  is  in  [0,  s  ],  p  (x)  is  the 
A 

a  a  a 

pressure  in  the  a  gear  gas  chamber  at  the  a  gear  stroke  x. 


p^  is  the  simple  graph  such  that  if  t  0,  p^  (t)  is  the  pressure  in  the  a 
^a  a 

gear  fluid  chamber  at  the  time  t. 
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P„  is  -  p,  [3^] 
a  a  a 


p_  is  the  simple  graph  such  that  if  t  0,  p„  (t)  is  the  pressure 
a  a 

in  the  a  gear  snubbing  chamber  at  the  time  t. 


Pj^  is  p^  -  p^  [s^]  for  the  type  I  strut  configuration  shown  in 

”3  a 

a  a 

Figure  B-1 . 


Pj^  is  p^  ■  strut  configuration  shown  in  Figure  B-3. 


BFU  ^  is  the  simple  graph  such  that  if  s  is  in  [0,  s„  ],  BFU  *^(s  )  is  the 

Si  3  w  &  21 

2 

shock  strut  upper  bearing  a  gear  reference  system  force  coefficient  for  a 
unit  bth  force  component  on  the  point  of  the  a  gear. 


BFL  ^  is  the  simple  graph  such  that  if  s  is  in  [0,  s„  ],  BFL  ^(s  )  is  the  a 

2  di  ^  Si  SL 

a 

gear  reference  system  lower  bearing  force  coefficient  for  a  unit  bth  force 

component  on  the  point  A  of  the  a  gear. 

a 


BMU  ^  is  the  simple  graph  such  that  if  s  is  in  [0,  s„  ],  BMU  ^(s  )  is  the  a 
a  a  b  a  a 

a 

gear  reference  system  upper  bearing  force  coefficient  for  a  unit  bth  moment 
component  on  the  axle  of  the  a  gear. 


BML  is  the  simple  graph  such  that  if  s  is  in  [0,  s  ]  BML  °(s  )  is  the  a 
<1  a  a  a 

gea'’  refer«*''ce  system  lower  bearing  force  coefficient  for  a  unit  bth  moment 


component  in  the  a  gear  reference  system  on  the  axle  of  the  a  gear. 
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If  the  a  gear  stroke,  s  ,  is  in  [0,  3„  ]  and  if  t  >  0,  the  shock  strut 

a 

force  ,  Fg  (t),  may  be  expressed  in  terms  of  the  forces  and  pressures  acting 
a 

on  and  in  the  shock  strut.  This  may  be  accomplished  by  isolating  the  piston 
of  the  type  I  shock  strut  as  seen  in  Figure  B-2.  From  this  figure  the 
following  expression  for  the  number  F  (t;  may  be  derived  from  a  summation 

O 

a 

of  the  forces  on  the  piston: 


=  PaMB  %  -  "t  -  “p 

a  a  a  a  a 

+  p  (t)  (A  -  A  (p  (t))) 

a  a  a  a 

-  Pa  <*A  "  V  -  V  ‘Ph,  > 

a  a  S  S  S  a 

a  a  a 


a  a  a 


With 

P.J,  (t)  Pj^  (t)  +  p^  (s^(t)) 
a  a  a 

and 

p  (t)  =  p  (t)  +  p  (s  (t)) 

^ 

a  a 

it  follows  that 


^'s  <P>  =  PaMB  *A  -  Ph  'P>  <''h  -  *P 
a  a  a  a  a 

*  P„  <t)  <V  -  V  'Ph„  '*•>>’  -  Pa  <Pa<‘»  *A  *  ‘'f 


where 


Ap  (t)  =  Aq  -  (s^(t)). 

a  a  a 
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When  the  shock  strut  has  bottomed  (i.e.  the  shock  strut  stroke  Is  less 
than  zero)  it  is  assumed  that  the  bottoming  load  is  derived  from  a  linear 
spring  and  a  linear  (viscous)  damper.  Thus,  for  this  case 

■  ^AMB 

a  a  a  a  a 

.  p  (t)  (A  -  A  (p  (t)))  -  Pj  (s^(t))  A,  .  Fp  (t) 

S  S  o  o  3  di 

a  a  a  a 

-'^0 

a  a 

It  is  now  required  to  derive  the  equation  for  the  pressure  in  the  shock 
strut  fluid  chamber  It  is  supposed  that  at  the  time  t  this  pressure  is 
directly  related  to  the  volumetric  expansion  of  the  hydraulic  fluid  and 
inversely  proportional  to  the  volume  in  the  fluid  chamber. 

This  relationship  may  be  expressed  by 

p  (t)  =  -  ig  (t)  if  g  (t)  is  <  0 

^a  (t) 

a 

and 

(t)  =  0  if  /Sg(t)  is  >  0. 

The  volumetric  expansion  at  the  time  t,  H  (t),  may  be  obtained  from  an 

'  a 

integration  of  the  sum  of  the  functions  that  define  the  fluid  volume  leaving 
the  fluid  chamber  through  the  orifice  and  the  rate  of  of  change  of  the 
volume  in  the  fluid  chamber.  Thus, 
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/6a'  =  S  *N  ''h  *  ’'h 


j  a 


a  a 


is  the  function  to  be  integrated  where  at  the  time  t  the  velocity  of  fluid 

through  the  main  orifice,  v  (t),  is  derived  from  the  principle  of 

H 

a 

conservation  of  energy  and  the  number  v  (t)  is  therefore 

h 

a 


Pu  (t) 


fcr  (t)  .  0, 
a 


V„  (t)  =  n 
ri 

a 


f  =  i  it)  0, 


Po  (t) 

V.,  (t)  =  -  (-  2  — # - 


for  Pj,  (t)  ■■  0, 


and  the  rate  of  change  of  the  volume  in  the  fluid  'unamber  is 


-  (A  -  A  (s  (t)))  3  (t)' 
H  P  a  a 

a  a 


A  similar  approach  is  used  to  determine  th';  pressure  at  the  time  t  in 
tr.e  snubbing  chamber.  The  requi-^ed  equations  ire 


P-T.  (t)  = 


/6s 


V„  (t)  a 

“s 

a 


P~  ( t )  r  0 

s 

a 


it  (t)  is  <;  0, 


if  /Sc,  (t)  >  0, 


where 


6s  ‘  V  %„  ‘fH.  V  *  V 


s  ”s  s 

a  a  a 
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Ph  (t) 
a 


v„  (t)  =  (2 


(t)  =  0 


-) 


0.5 


Pu  (t) 
a 


v„  (t)  =  -  (-  2 


>0.5 


for  p  (t)  >  0, 
a 


for  p  (t)  =  0, 
a 


for  p  (t)  <  0, 

h 

a 


V„  (t)'  =  A„  s^'(t). 
a  a 


It  is  assumed  that  the  gas  volume  in  the  shock  strut  at  the  time  t  is 
defined  by 


’a  =  V  -  % 

a  0  a 

a 


It  is  further  assumed  that  n  is  a  positive  number  such  that  the 
thermodynamic  process  of  the  gas  in  the  shock  strut  gas  chamber  during 
compression  is  represented  by 


a 


Pa  =  Pa  ^ - 


For  the  type  II  shock  strut  at  the  time  t  shown  in  Figure  B-3  the 

pressures  p  (t),  p  (t)  and  p.  (t)  are  calculated  from  the  relations 
a  »a 

derived  for  the  type  I  shock  strut.  There  is,  however,  a  different 

equilibrium  equation  from  which  the  number  F^  (t)  is  derived.  The  following 

a 


99 


expression  for  F  (t)  may  be  derived  from  Figure  B-4: 
a 


"s  '  “amb  *a  -  »T  “b  -  *P  V 

a  a  a  a  a  S  S 

a  a 

+  p  (t)  (A  -  A  (p  (t)))  -  p  (s  (t))  A  (s  (t))  +  F  (t)  for  s  (t)  >  0, 

”s  *^3  "S  a  ^  ^a  “  *^a  ^  ~ 

a  a  a  a 

and  since 

Pt  (t)  =  p^  (t)  +  p^  (s^(t)), 

^a  a  a 

Pt  (t)  =  Py  (t)  +  p^  (t), 

S  S  a 

a  a 

and 


A  -  A  -  A 
A  H  H_  ’ 
a  a  3 

a 

it  follows  that 


^3  “  °AMB  *A  “  ^ 

a  a  a  a  a  S 

a 


*  ph,  <V  -  \  'X  **  * 

3  3  3  3  a  a  a 

a  a  a  a 


for  s  (t)  >  0. 
a  — 


For  the  case  where  s  is  less  than  zero,  the  same  stiffness  and  damping 

a 


terms  that  were  used  for  the  type  I  shock  strut  are  to  be  included. 


Also,  the  rate  of  change  at  the  time  t  of  the  volume  in  the  fluid 
chamber  is 


Vh  (t)' 


(t)) 


(t)  -  (Aj^ 

a 


A  (s  (t)))  s  '(t). 
r  a  a 

a 


To  compute  the  shock  strut  friction  force  function  it  is  first 
neccessary  to  calculate  the  force  and  moment  functions  for  the  gear  axle. 
This  may  be  done  using  Newtonian  mechanics.  From  this  method  it  seen  from 
equation  (C-1)  that  the  axle  force  for  the  a  gear  may  be  expressed  by 


ICX) 


"  JVu  ;I  ^;I  ^G  ' 

A  J  M  a’t’ta  a  a 

a  a 


The  a  gear  axle  velocity  function  in  the  parenthesis  of  this  equation 
may  be  found  in  an  expanded  form  in  equation  (7)  in  the  Derivation  of  the 
Equations  of  Motion  section  in  the  main  body  of  the  report. 


If  the  indicated  differentiations  and  integrations  are  performed  amd 
the  terms  involving  the  rates  of  change  of  gear  and  airframe  deformations 
are  omitted,  then  the  following  expression  for  the  numbers 


F  (t);  b  =  1,  3  is  obtained: 

% 


Fq  (t)  =  -  Qq  \t)  (Vg^'(t)  +n/(t)  Vg3(t)  Vg^(t)) 


a  a 


“  ”w  <^G  +ne^(t)  Vg\t)  Vg3(t)) 


a  a. 


-  Qq  (Vg'^'(t)  ♦ng’(t)  Vg'^(t)  -Ug^Ct)  Vg'(t)) 


a  a. 


-  "w  Qg  ®<t)  (Ij  ^ 


-  gj  '■  Qg  ’(t)  *  Og  '’(t))  Sg  Sgh 

a  a^  a  a^. 


Qq  Sp^Hg^ct)  Hg^ct)  (1^  d 


a  a. 


*  g,  “Qg  "(t)  *  ’(3,(1))  00"  (t))  §'•3 

H  3  3  3 

U  V 


e 

a 


Hg  ’(t)  (g,  ^  '’(s^(t))  Sa  5^^ 

a  a  a  “ 
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-M,  (g/ . 


bs 


2M„  a,«(t)  S,p  '■•(s.tt))  3^'(t)  Qo  "(t)  8,,  8,, 

a  a^  a  a^ 


(t)  f.  ‘>•(3,(1))  3,’(t)  8,,  8^^ 


W  G  A  a  a 

a  a  a 


h„ 

J  A 


\  =A  '"a 
a  a 


(3.(t))  (s.'(t))‘  8_  -  M„  i.  '■'( 


bh-'Sl  5,  •■■(3^(t))  3^-(t)  8,, 
a  a 


-  'Sf  ‘t’BG.  <9g  Ib'’"*'’  "(t)  8 


a  A 


bh 


-  (t))  ‘^"(t)  8 


a  A 


bf 


r,  «(t)  7  =(t)  a,,  “(t)  8,,  *  «„  g  7b!'^)  Qg 


The  contribution  of  the  axle  moment  from  the  acceleration  of  the  wheel 
mass  typically  has  a  negligible  influence  on  the  strut  friction  forces. 
Therefore,  this  moment  contribution  will  not  be  formulated.  However,  the 
axle  moment  from  the  ground  forces  may  in  some  cases  have  a  significant 
effect  on  the  strut  friction  force.  It  is  assumed  that  the  axle  moment  from 
the  ground  force  is  parallel  to  (t). 


Thus,  the  axle  moment  from  the  ground  force  is 


(iy  (t)  •  (w^  (t)  +  v„  (t))  X  F,  (t))  i,,  (t) 


W 


Therefore,  It.  is  f'cind  that 
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Mq  (t)  =  ^t)  (w^  °(t)  >  Vy  (t))  Fq  ^t) 


The  moment  components  at  the  time  t  in  the  a  gear  reference  system  are 
determined  from 


Mq  (t)  =  Mq  (t)  (t)  *  Jq  (t) 

a.  a. 

s  ^ 


(t)  =  (t)  “(t)  y„  '(t) 


s  °  “d 


The  next  step  in  the  derivation  of  the  shock  strut  friction  force  is  to 

determine  the  bearing  reactions  between  the  piston  and  the  cylinder.  When 

computing  these  reactions  the  assumption  will  be  made  that  the  forces  on  the 

piston  (and  axle)  from  the  acceleration  of  this  mass  can  be  omitted.  This 

permits  the  bearing  reactions  to  be  computed  from  a  simple  expression 

involving  the  axle  forces  and  moments.  Therefore,  it  follows  that  the 

number  F  (t)  can  be  determined  from  the  following  equation: 
r 

a 


3  '  (t) 

^  b  c 

F„  (t)  - - ,  (Mnn  (BFU^°(3  (t))  F_  (t)  +  BMU^^(3^(t))  M„  (t)) 

^b  “c 


-»■  U_,  (BFL  (3^(t))  F„  (t)  +  BML^®(s  (t))  M„  (t))) 
'  BL  a  a  G .  a  a  G . 
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p.  (3  (t)  A  (s  (t))  (t) 

A  a  Pa  o 
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APPENDIX  C 


MODIFICATION  OF  LAGRANGE'S  EQUATIONS  FOR  NON-HOLONOMIC  CONSTRAINTS 

AND  QUASI-COORDINATES 

3 

Suppose  that  V  is  a  volume  in  E  space  containing  a  collection  of  mass 

particles  and  that  S  is  the  surface  of  this  volume.  Further,  suppose  that  Q 

is  a  point  fixed  in  space  or  moving  with  a  constant  velocity  and  that  r  is  a 

12  3 

vector  function  such  that  if  t  0  and  if  (x  ,  x  ,  x  )  is  a  point 
-  1  2  3 

in  V,  r(x  ,  x  ,  x  ,  t)  is  the  vector  at  the  time  t  from  Q  to  the  point 
labeled  (x\  x^,  x^)  in  V. 

Now  suppose  that  5  r  is  the  vector  function  such  that  if  t  0  and  if 

123  Sv-123 

(x  ,  X  ,  X  )  is  a  point  in  V,  Or(x  ,  x  ,  x  ,  t)  is  a  vector  which  at  the 

time  t  is  tangent  to  the  constraints  (which  are  assumed  to  be  fixed  at  the 

12  3 

time  t)  on  the  point  labeled  (x  ,  x  ,  x  )  in  V  but  otherwise  arbitrary  in 
magnitude  auid  direction. 

Further  suppose  that 

-  12  3 

P  is  the  vector  function  such  that  if  t  0  and  if  (x  ,  x  ,  x  )  is  a  point 

in  V,  P(x  ,  X  ,  X  ,  t)  is  the  body  force  per  unit  volume  at  the  time  t  at 

the  point  labeled  (x\  x^,  x^)  in  V, 

“  12^ 

Z/ is  the  vector  function  such  that  if  (x  ,  x  ,  x'^)  is  a  point  on  the  surface 
—  1  2  3 

of  V,  l/(x  ,  X  ,  X  )  is  the  unit  vector  normal  to  the  surface  at  the  point 
12  3 

labeled  (x  ,  x  ,  x  )  on  the  surface  of  V. 
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is  the  vector  function  such  that  if 
_  1  2  3 

on  the  surface  ofV,T  (x,x,x,t) 
at  the  time  t  on  the  surface  point  (x\ 
x^  x^), 


1  2  3 

t  0  and  if  (x  ,  X  ,  x  )  is  a  point 

is  the  surface  force  per  unit  area 
2  3 

x  ,  X  )  where  the  normal  vector  is 


U  123  1  2  3 

I  is  the  simple  surface  sucn  that  if  (x  ,  x  ,  x  )  is  in  V,  '  (x  ,  x  ,  x  ) 

is  the  mass  density  of  the  collection  of  mass  particles  at  the  point  labeled 

(x\  x^,  x^)  in  V,  and 


Sw  is  the  simple  graph  such  that  if  t  0,  Sw(t)  is  the  (virtual)  work 
done  at  the  time  t  by  the  component  of  the  body  and  surface  forces 
in  the  direction  of  the  vector  Sr(t)  acting  through  the  distance  lSr(t)|  , 
then  from  Reference  ( 1 ) ,  the  Principle  of  Virtual  Work  states  that 


(C-1) 


The  principle  of  Virtual  Work  may  be  used  as  a  basis  for  the  derivation 
of  Lagrange's  equations.  However,  before  Lagrange's  equations  can  be 
derived,  the  vector  function  Sr  must  be  rewritten  in  terms  of  the 
coordinates  (i.e.  functions)  that  describe  the  motion  of  the  collection  of 
mass  particles. 


A  vector  function  that  could  be  used  for  this  purpose  is  a  modification 
of  the  velocity  vector  function.  To  develop  this  vector  function  suppose 
that  n  is  a  positive  integer,  each  of  q  ;  a  =  1 ,  n  is  generalized  coordinate 
and  R  is  a  vector  function  of  class  C  such  that  if  t  >  0  and  (x  ,  x  ,  x^) 
is  in  V, 
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r(x\  x^,  x^,  t)  =  S(x\  x^,  x^,  t,  q®(t);  a  =  1,  n) 
at  the  time  t. 


The  velocity  vector  is  therefore 


r  (x\  x^,  x^,  t)  =  R  b(x\  x^,  x^,  t,  q*(t);  a  =  1,  n)  q^’(t) 


+  R  (x\  x^,  x^,  t,  q^(t);  a  =  1,  n). 
’■^t 


In  this  expression  the  R^^  term  has  a  contribution  to  the  velocity  for 

’^t 

the  case  where  the  constraints  are  time  dependent.  However,  for  the 
calculation  of  Sr,  this  contribution  is  not  included  since  the  constraints 
are  fixed  when  S r  is  determined. 


Therefore,  the  expression 

Sr(x\  x^,  x^,  t)  =  R  b(x\  x^,  x^,  t,  q*(t);  a  =  1 ,  n)  q^'(t)  (C-2) 

I  H 

may  be  used  to  derive  the  vector  function  Sr  that  could  be  used  in  equation 
(C-1)  where  it  is  understood  that  (within  the  constraints)  that  each  of  the 
coordinate  derivatives  are  independent  and  arbitrary  in  magnitude. 

This  substitution  for  the  Sr  in  equation  (C-1)  will  provide  the 
formulation  needed  for  the  derivation  of  Lagrange's  equations  for  the  case 
where  the  motion  is  expressible  in  terms  generalized  coordinates. 

For  this  purpose  suppose  that  T  is  a  simple  graph  such  that  if  t  >  0, 
T(t)  is  the  kinetic  energy  of  the  collection  of  mass  particles  at  the  time 
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t.  Further,  suppose  that  K  is  a  simple  surface  such  that 

A 

T  =  KEq'^;  c  =  1,  n,  q°';  c  =  1,  n] .  (C-3) 

Therefore,  if  each  of  Q^;  a  =  1,  n  is  a  simple  graph  such  if  t  0, 

Q  (t);  a  =  1,  n  are  the  generalized  forces  acting  on  the  collection  of  mass 

ci 

particles  at  the  time  t  and  the  acceleration  function  term  is  expressed  in 
terms  of  the  kinetic  energy  function  as  illustrated  in  Reference  (5),  then 
the  Principle  of  Virtual  Work  may  be  rewritten  in  the  form  of  power  rather 
than  work  as  follows: 

((K.^a,)’  -  K.^a  -  Q^)  q*’  =  0.  (C-4) 

Lagrange's  equations  may  be  easily  derived  from  equation  (C-4)  for  the 
case  where  the  velocity  functions  are  independent. 

In  some  cases  it  is  preferable  to  write  the  equations  of  motion  in 
terms  of  non-general ized  (quasi)  coordinate  velocity  component  functions. 

For  example,  it  is  frequently  desirable  to  have  the  rigid  body  motion  of  a 
collection  of  mass  points  be  determined  from  the  body  axis  components  of  the 
velocity  vector  function  and  the  angular  velocity  vector  function.  These 
are  not,  however,  generalized  coordinate  velocity  component  functions  since 
the  position  of  the  body  cannot,  in  the  general  case,  be  determined  from  an 
integration  of  these  velocity  component  functions. 

Suppose  that  each  of  V  ;  a  =  1 ,  n  is  a  quasi-coordinate  velocity 
component  function  and  that  each  of  A^;  a  =  1 ,  n;  b  =  1 ,  n  is  a  simple 
surface  of  class  at  each  of  its  points  and 
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(C-5) 


..3.  ,.3rC  ^ 

V  =  k^iq  ;  c  =  1 ,  n]  q  ' . 


Further,  suppose  that  each  of  Bj^;  a  =  1 ,  n;  b  =  1 ,  n  is  a  simple 
surface  such  that 


q^'  =  B^[q°;  c  =  1 ,  n] 


It  follows  then  that  if 


a 

and 


Ajj[q  ;  c 


=  1 ,  n] 


^  =  Bj^[q  ;  c  =  1,  nj, 


then  is  unity  for  the  case  where  a  is  equal  to  c  and  zero  for  the 

case  where  a  is  unequal  to  c.  Also,  it  is  seen  that 


'■  =/S^' 


(C-6) 


If  equation  (C-6)  is  substituted  into  equation  (C-4)  and  if  each  of 
S^;  b  =  1,  n  is  a  simple  graph  such  that 


\  =  /Sb  “b- 


(C-7) 


then 


</Sb 


V  '' 


=  0. 


(C-8) 


It  remains  now  to  replace  K_  a,  and  K,  a  with  functions  expressible  in 

J  Q  *  Q 


113 


terms  of  quasi-coordinate  velocity  component  functions.  For  the  purpose, 


suppose  that  K  is  a  simple  surface  such  that 
T  =  K[q°;  c  =  1,  n,  V°;  c  =  1,  n] . 


(C-9) 


The  kinetic  energy  functions  in  equations  (C-3)  and  (C-9)  may  be 
differentiated  to  yield 


T'  =  K  a  q^'  +  K  a,  q*" 

;q  ;q  ’ 

=  K  a  q^*  +  K.„a  V^'. 

f  q  >  • 


(C-10) 


But,  from  equation  (C-5) 


-.a .  .a  b.c,  n a^  d  ,  ,  i  b_ 

V  •  =  Ac  q  ’  q  '  +  A  [q  ;  d  «  1,  n]  q 
D ,  q  D 


Therefore,  by  comparison  of  like  term®  in  equation  (C-10)  it  is  seen 


K.^a,  =  K.yd  A^[q°;  c  =  1,  n]. 


K  a  =  K  a  +  K  „d  aJ  a  q^' 

;q  ;q  ;v  b;q 


If  these  expressions  are  substituted  into  equation  (C-8)  then  it  is 


found  that 


^  (K  ,,d)'  +  K  „d  V* 

A-'b  a  ;V  ;V  a;q  A-' e 


K„a  -  K  .,d  a'^  a  V*]  =  s  . 

>  q  ;  V  c ;  q  e  ■'  b  ’ 


and  if 
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-  *a,q=>/Se 


then 


((K.,b)'  -nj  K.,d  K.^a  -  Sj,)  v”  = 


(C-11) 


(C-12) 


is  the  Principle  of  Virtual  Work  expressed  in  terms  of  quasi-coordinate 
velocity  component  functions. 


Now  suppose  that  m  is  a  positive  integer,  but  less  than  n  and  that  each 
of  a  =  1 ,  m  is  a  simple  surface  such  that  if  t  0, 

r^bCt,  q°(t);  c  =  1,  n)  V^’Ct)  =  0  (C-13) 

is  a  constraint  on  the  collection  of  particles  at  the  time  t.  Note  that  the 
constraint  relation  is  linear  in  terms  of  the  quasi-coordinate  velocity 
component  functions. 

Further  suppose  that  F  is  a  simple  surface  such  that  if  t  0, 

^  . 

F  (t,  q°(t);  c  =  1,  n,  V  (t);  b  =  1,  n)  (C-14) 

di 

at  the  time  t  and  if  each  of  ^  n  is  a  simple  graph  (a  Lagrangian 

multiplier)  such  that 

X^(t)  F  (t,  q°(t);  c  =  1,  n,  V^(t);  b  =  1,  n)  =  0  (C-15) 

3 

at  the  time  t. 


It  follows  then  that  for  the  case  of  non-holonomic  constraints 
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expressed  in  the  form  of  equation  (C-13),  Lagrange's  equations  may  be 
written  as  follows: 


-n' 


K.v<l 


K  a 
b  ;q 


"b  " 


(C-16) 


The  use  of  equation  (C-7)  to  find  the  S^;  b  =  1,  n  functions  is  often 

difficult.  A  more  direct  approach  is  to  develop  another  candidate  for  Sr 

in  equation  (C-1)  in  terms  of  the  quasi-coordinate  velocity  component 

functions.  Suppose  that  A is  a  vector  function  such  that  if  t  >  0  and 

b  — 

1  2  3 

(x  ,  X  ,  X  )  is  in  V,  then  the  the  sum 

Aj^(x\  x^,  x^,  t,  q^(t);  a  =  1,  n)  V^(t)  (C-17) 

S-  1  2  3 

r(x  ,  X  ,  X  ,  t)  at  the  time  t  for  the  point  labeled 

12  3 

(x  ,  X  ,  X  )  in  V.  Note  that  the  vector  in  (C-17)  is  the  same  as  the  vector 
in  (C-2).  However,  the  vector  sum  in  (C-17)  arises  naturally  from  the 
differentiation  of  the  vector  function  r  when  quasi-coordinate  velocity 
component  functions  are  used.  It  follows  then  from  equation  (C-1)  that 


A  )  dV  + 

b 


(f^  '  A.  )  dS. 

b 


(C-18) 
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APPENDIX  D 


DEFINITION  OF  THE  MASS  TERMS 

The  mass  terms  which  appear  in  equation  (8)  are  defined  below.  These 
terms  are  derived  from  the  operations  required  when  equations  (5),  (6)  and 
(7)  are  substituted  into  equation  (1). 


(It  ®  *  a„  '  (u/  *  H/fl,  1,  I.  2,  I,  3,  =J)) 

®  a  a 


=a.  ^  ^  *0  ’‘g  “g  ^ 

f  a  a  a 


<1g'’‘«a‘’*  *■«; 

a„  a  a  a 

f  m 


')) 


“  a_  a  a  a  <=>  „  , 

I  m 


a  a 
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(M  ^ 


3  3 

a  a 

(M  )[3 J 
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\  a  ®  a  ^  ^ 

a  b 


X  ( 
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O 


■A  ‘-a'  •  “g  ''“al  "a  >  0,,  dV„ 

a^  a  a 
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% 

b  Oq  >  ^ 

a 

=  /v„4>b>b'  5,e  P, 


B 


N, 


G 

>  I 


/, 
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<«  °'b  \)t=a'  6a  ' 

a 
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Vg  % 

b  0G  ;  a  =  1.  Nq] 
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B  c 
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ad  0  a 


f,  dV^ 

G  u 
a  a 


”g  r 
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^B'^W 

(M  ^  ^[s^,  ] 
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b  Oq  >  ^  Nq] 

a 

=  /v„\b3  4''‘t’B‘'Sl?  Pb'”! 
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a 
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a  a  G  “G  ''G 
n  c  a  a  a 


G  '^^G 
a  a  a 


J  V  ®hb3  ^G  ^®A  ^A  *-®a^  ^G 

a  a  a  a  a 

am  g 


X  a,,  '■(fco  "[3^,  Po  ]  5?  P„  <!»„ 


a  A 
n  a 


a  a 


])) 


123 


(M  ®  ^)[3  J 
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'^G  *^0 
a  a 


(M 
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b>IS-  I 
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APPENDIX  E 


DERIVATION  OF  THE  GROUND  FORCE  EQUATIONS 

The  derivation  for  the  component  of  the  ground  force  on  the  tire  that 
is  normal  to  the  ground  plane  was  given  in  the  body  of  the  report.  It  was 
established  as  being  functional  dependent  only  on  the  simple  graph  v  .  The 

n 

a 

ground  force  components  in  the  ground  plane  are  somewhat  more  complicated. 
One  of  the  complications  is  that  the  ground  forces  are  dependent  on  whether 
or  not  the  tire  footprint  is  slipping  on  the  ground  surface.  Prior  to  "spin 
up"  of  the  wheel  the  tire  footprint  is  slipping  on  the  ground  surface  and 
the  force  in  the  ground  plane  is  in  the  direction  opposite  to  the  velocity 
vector  of  the  tire  footprint.  The  magnitude  of  this  force  is  determined  by 
multiplying  the  magnitude  of  the  ground  force  component  normal  to  the  ground 
plane  by  the  friction  coefficient  which  may  be  defined  in  the  following 
manner . 

Suppose  first  that  R  is  the  simple  graph  such  that  if  t  >  0,  R  (t)  is 

3l  —  3 

the  slip  ratio  of  the  a  gear  tire  at  the  time  t  (i.e.  the  ratio  of  the 

magnitude  of  the  tire  footprint  skidding  velocity  to  the  magnitude  of  the 

velocity  of  the  axle  parallel  to  the  ground).  Also  suppose  that  fJi^  is  the 

a 

simple  graph  such  that  if  x  ^  0,  (x)  is  the  friction  coefficient  between 

a 

the  a  gear  tire  and  the  ground  at  the  slip  ratio  x. 

After  spin  up  (i.e.  the  slip  ratio  is  zero)  the  forces  on  the  tire  in 
the  ground  plane  change  their  functional  dependence.  The  drag  force 
component  on  the  tire  in  the  ground  plane  is  assumed  to  be  a  number  times 
the  normal  ground  force  component.  The  side  (or  cornering)  force  component 
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is  defined  through  the  following  functions. 


Suppose  that  Is  the  simple  graph  such  that  If  t  >0,  (t)  Is 

a  a 

the  cornering  angle  of  the  a  gear  tire  at  the  time  t.  Suppose  further  that 

Is  a  simple  surface  such  that  the  point  (x,  y,  (x,y))  belongs  to  F^ 
a  a  a 

only  If  at  a  cornering  angle  x  and  a  normal  ground  force  component  y, 

F^  (x,  y)  Is  the  side  force  component  on  the  tire, 
a 


Also,  at  a  time  t  after  spin  up  of  the  a  gear  wheel  the  component  of 
the  ground  force  parallel  to  the  unit  vector  1^  (t)  (l.e.  the  drag  force 

component  on  the  tire)  Is  determined  by  multiplying  the  normal  ground  force 

component  by  the  rolling  resistance  jJL^  which  Is  assumed  to  be  a  number. 

a 


The  vector  functions  needed  to  define  the  kinematic  relations  are 
described  as  follows: 


r^  *  Is  the  vector  function  such  that  If  t  0,  r^  '(t)  Is  the  velocity  of 
a  a 

the  a  gear  axle  point  A  with  respect  to  the  point  Q  at  the  time  t. 

3 


'  Is  the  vector  function  such  that  If  t  0,  rp^  '(t)  Is  the  velocity  of 
a  a 

the  a  gear  tire  footprint  reference  point  C  with  respect  to  the  point  A 

F  a 

a 

at  the  time  t. 


Each  of  ly  ;  b  = 
% 

ly  (t)  is  1^, 

^3 


and 


1,  3  is  a  unit  vector  function  such  that  if  t  0, 

(t)  is  the  cross  product  of  k^  (t)  (see  Appendix  G)  and 

*2 

(t)  is  (t)  X  (t)  at  the  time  t. 
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The  derivation  of  the  kinematic  relations  may  now  be  accomplished  with 
the  help  of  the  unit  vector  systems  and  the  transformations  defined  in  the 
main  body  of  the  report  and  Appendix  G. 

The  a  gear  right  hand  orthogonal  unit  vector  functions  ;  b  =  1 ,  3 

^b 

for  the  wheel-ground  reference  system  may  be  expressed  by  the  following 
functional  relations: 


X  I 


3 


f_  X  (k„  X  i.) 


|f  X  (ic^  X  I  )| 

^2 

/So  ‘’Qo  °  (7b'  i,  *  7b^  Ij) 
[(^0  "Qb  °  7b')"  -  (^G  "aG  ° 


=  I 


3' 
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The  axle  velocity  of  the  a  gear  may  be  obtained  from  the  terms  in 
equation  (7)  in  the  main  body  of  the  report.  This  velocity  is  expressed  by 


+^2Q(t)  X  (I^  (t)  +  (t)  +  (s^(t),  t)) 

a  a  a  a 


*  I 


+  0^2  (t)  X  (g^  (t)  +  (3^(t),  t)) 

a  a  *a 


-  I. 


On  ^4  •«,  (s  (t),  t)  3„'(t) 


A  ;s  a 
a’  a 


BG»  ■  a 

A  a 

^b 

+  <i)G  (t),  t)  Qq  ^'(t), 

A.  a  a 


or  in  expanded  form 


a  f  f  a 

.  1  (t)  /‘»'n/(t)  (t)  (g,  = .  =(A^(t))  8.P  S,, 

f  a  a  a 

c 

»  (t)  (t)  (gj  °  .  I,  =(3^(t))  Spr 

a ,  a  a  a 

d 
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a 


+  3„  (t)ct)  "^(3  (t),  (t))  q_  ‘’-(t). 


'G  '*"G,  '  a 

A 


Thus ,  with 


G  '  "  ^G 
a  a 


(t)  l„  (t), 


‘b  ^ 


the  axle  velocity  is 
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'(t)  =  ly  (t)  Jy  °(t)  (Vg^(t)  +  e^P^n3®(t)  (1^ 

a  a  f  a.  a 

c  b 

*  “(t)  (g,  o . 

a  a  a 

q 

*  Qg  "(t)  '^BJt)  T„  <‘>  'Sa  ”  *  S,f 

3l  ,  D  3^3 

d  h 

*  To  Te^t)  T„  “(U'J^BG  <‘»  ^b'’''^* 

a  b  a  a  a 

d  h  f 

^G  (t)  (t)  *^W  (t) 

b  a.  a 

a  n 

(t)  (t)  (t)*i’G  &G  '’g 

a^  b  a,  A  a  a 

d  ha. 

b 

With  the  assumption  that  the  stroking  velocity  has  a  negligible  effect 

on  the  velocity  r_.  '(t),  this  vector  may  be  approximated  by 
r  A 

a 

TpA  '(t)  =  (iy  (t)  (r2g(t)  X  (t)  +  (HgCt)  ■  kg  (t))  icg  (t)  X  Vy  (t) 

a  a.|  a  a 

+  (t)  X  Wg  (t)  +^7^  (t)  X  Wg  (t)  +i7y  (t)  X  Vy  (t))  '  ly  (t)) 

a  a  a  a  a  a 

+  (ly  (t)  (i7g(t)  X  Wg  (t)  +  (t)  X  Wg  (t))  ’  ly  (t)) 

A^  a  a  a  a^ 

+  iy  (t)  Vy  '(t). 

^3  ^ 


After  the  indicated  vector  operations  are  performed  it  is  found  that 

'(t)  =  i„  (t)  an^\t)  “(t)  /e”  y/(t)  s,,  Tw  '(t)  8p 

a  a,  a  b  *  a 

1  r 

+  yS  '"(t)  a.  “^(t)  ShH^ 


(Vy  (t)  ySc  ^(t)  Qg  ®(t)  ^g^t)  '(,)  g 
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(t)  °(t)  Qq  ‘^(t)  y/ 

a  a-  d 


,  '“(t)  y^ht) 

’^3  " 


8.8, 


+  (Hy  (t)  °(t)  ‘'(t)  Qq 

a  a  a^ 
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a  a  a^  a^ 


‘'(t)  Qq  "^(t)  J/(t) 
%  d 


y„  ^t)  S«) 

a 

Sfs  S,t  T„  '(t)  S® 

a 
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*  '^w  (t)  (t)  ‘'(t)  Qq 

a  a  a^ 

+  (t)  (cr^g^’ct)  w^  °(t)  y^ 

"a„  a  1 


(t)  (cr^g^’ct)  w^  °(t)  yB'^(t) 
a^  a  b 

(t)  w^  °(t)  a.  ^(t)  yg^ct) 
^a  a  a 2  d 


+  !y  (t)  Vy  ’(t). 

as 


b  g 

(t)  y^^t)  e'’=3  y^  '(t)  S®)) 

d  a 

g 

,V.  Z®*'  S,,S,,T„^(t)8?) 

D  a 


a 

m 

Sfs  8=t  )'«'<«  8^)) 

a 

m 


The  a  gear  wheel  has  reached  a  spin  up  condition  at  the  time  t^^  if 


"■a  '<‘su>  *  '■fa  '‘'su»  ■  '‘su>  '  “• 

a  a  a^ 


For  t  <  t  the  ground  force  parallel  to  the  ground  plane  is  in  the 
oU 

direction  of  the  vector 


(-  r^  *(t)  •  ly  (t)  -  '(t)  '  (t))  (t) 

a  a^  a  a^  a^ 

+  (-  *(t)  •  (t)  -  '(t)  •  (t))  (t) 

a  a^  a  a^  a^ 


A(t)  ly  (t)  -  B(t)  (t). 
^  ^2 


Now  suppose  that  I/'  is  the  simple  graph  such  if  t  is  in  [0,  t  ), 

a  o  u 


I/'^Ct)  =  tan 


-1  B(t) 
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at  the  time  t. 


Therefore,  since  the  slip  ratio  at  the  time  t  is  defined  by 
(?.  ’(t)  +  ’(t))  ■  L 


FA  W 

a  a 


R  (t)  = 
a 


1 


a  a 


and  since  the  ground  force  vector  may  be  expressed  in  terms  of  its 
components  in  the  wheel  plane  reference  system  by 


a  a  a. 

b 


it  follows  that 


’  =  -  Mo  '"a’  ^0  ^ 

a  a  a 


^G  "  =  -  Mo  ^G  ^ 


''c  ^  =  "g 

a  a 


For  t  >  tgy  the  ground  force  component  parallel  to  (t) 

is  derived  from  the  wheel  and  the  tire  rolling  resistance.  The  ground  force 
component  parallel  to  I,,  (t)  is  derived  from  the  tire  cornering  effect. 

If  each  of  C  and  D  is  a  simple  graph  such  that  if  t  >  t^^, 

c(t;  =  Mt)  •  i^,  u), 

a  a, 

and 
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D(t)  =  '(t)  *  (t) 

a 

at  the  time  t,  the  cornering  angle  of  the  a  gear  tire  at  a  time  t  may  be 
approximated  by 


(t) 


tan 


-1  D(t) 
C(t) 


It  follows,  therefore,  that 
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APPENDIX  F 


DERIVATION  OF  THE  ARRESTING  HOOK  FORCE  EQUATIONS 

The  following  derivation  of  the  arresting  hook  force  equations  will 
accornmodate  yawed,  rolled  and  drift  landings  in  addition  to  landings  where 
the  arresting  hook  engages  off  the  midpoint  of  the  cable.  In  this 
derivation  it  is  assumed  that  the  arresting  hook  and  the  cable  are  massless. 
Therefore,  the  orientation  of  the  arresting  hook  and  the  cable  may  be 
completely  determined  from  the  dynamics  of  the  airframe  and  the  geometry  of 
the  arresting  hook  and  the  cable  arresting  system.  The  arresting  force  from 
the  cable  is  assumed  to  be  a  function  of  the  cable  runout  only.  It  is 
further  assumed  that  there  is  friction  between  the  arresting  hook  and  the 
cable  interface.  The  derivation  of  the  equations  for  the  arresting  hook 
forces  accounts  for  both  slipping  and  not  slipping  of  the  arresting  hook  on 
the  cable. 

Some  of  the  equations  involved  are  nonlinear  algebraic  and 
trigonometric  relations  for  which  a  closed  form  solution  is  not  viable.  For 
these  cases  an  iterative  approach  is  described  that  may  be  used  to  obtain  a 
solution . 

The  points  on  the  arresting  hook,  cable  and  the  ground  that  need  to  be 
identified  ire  described  as  follows: 

H  is  the  point  of  attachment  of  the  arresting  hook  to  the  airframe. 

HP  is  in  the  arresting  hook  plane  of  symmetry  and  at  the  center  of  the  cable 
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cross  section  subsequent  to  cable  engagement. 

HC  is  the  point  on  the  ground  such  that  the  straight  line  perpendicular  to 
the  ground  containing  the  point  HP  contains  the  point  HC. 

LP  is  the  lateral  pivot  point  in  the  arresting  hook  shank. 

SL  is  the  point  of  exit  of  the  cable  from  the  left  hand  cable  sheave. 

SR  is  the  point  of  exit  of  the  cable  from  the  right  cable  sheave.  The 

straight  line  containing  SL  and  SR  is  assumed  to  be  parallel  to  the  straight 
line  containing  the  ground  fixed  unit  vector  which  is  defined  in  the  main 
body  of  the  report. 

S  is  the  midpoint  between  SL  and  SR. 


Q  is  the  ground  reference  point. 


N  is  the  airframe  reference  point. 


The  numbers  required  for  the  calculation  of  the  arresting  force  are 
defined  as  follows: 

t^  is  the  time  at  which  the  arresting  hook  engages  the  cable. 


'^BAR  horizontal  travel  (parallel  to  1^  of  the  point  HP  from  time 


equal  to  zero  to  time  equal  to  t^. 
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is  the  maximum  runout  of  the  arresting  cable  system. 


CD  u  is  initial  (i.e.  at  a  time  t  before  contact  with  the  ground  or 

r  Hq 

cable  engagement)  angle  between  the  arresting  hook  shank  and  airframe  unit 
vector  i_  (t). 


is  the  friction  coefficient  between  the  cable  and  the  arresting  hook. 


be  is  the  cable  semispan  (half  of  the  distance  between  the  points  SL  and 
SR). 


The  following  vector  is  required  to  establish  the  cable  geometry: 

Vg  is  the  vector  parallel  to  from  the  point  S  to  the  ground  and  Vg  is  the 

minimum  distance  from  the  point  S  to  the  ground.  Therefore, 


=  V, 


.3. 


The  following  vector  functions  are  required  in  the  derivation  of  the 
arresting  hook  force  equations: 

I  ,  I  and  I„  are  orthogonal  unit  vector  functions  such  that  if  t  >  0, 
I„  (t)  and  I„  (t)  are  fixed  in  the  plane  of  the  cable  with  I„  (t)  =  I„  , 

Hi  2 

I  (t)  is  normal  to  the  plane  of  the  cable  and  directed  down  relative  to 

”3 

the  pilot  and  I  (t)  =  L  (t)  x  L  (t)  at  the  time  t. 

1  “2  ”3 

1  is  the  vector  function  such  that  if  t  >  0,  1.  (t)  is  the  vector  from  the 
L  ■“  Lj 
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point  HP  to  the  point  SL  at  the  time  t. 

1  is  the  simple  ^raph  such  that  if  t  >  0,  1, (t)  is  the  magnitude 
of  the  vector  at  the  time  t. 

is  the  vector  function  such  that  if  t  >  0,  I_(t)  is  the  vector  from  the 
point  HP  to  the  point  SR  at  the  time  t. 

is  the  simple  graph  such  that  if  t  >  0,  l_(t)  is  the  magnitude 
of  the  vector  lp(t)  at  the  time  t. 

itO  is  the  vector  function  such  that  if  t  0,  M(t)  is  the  vector  parallel 
to  the  unit  vector  (t)  from  the  straight  line  containing  the 

points  SL  and  SR  to  the  point  HP  at  the  time  t. 


RO  is  the  simple  graph  such  that  if  t  0,  R0(t)  is  the  magnitude 
of  the  vector  R0(t)  at  the  time  t. 


R„  is  the  vector  function  such  that  if  t  >  0,  Ru(t)  is  the  vector  from  the 
M  —  H 

point  S  to  the  point  LP  at  the  time  t. 


is  the  vector  function  such  that  if  t  0, 
S„(t)  =  (R„(t)  •  I„  (t))  I„  (t). 


4 


is  the  simple  graph  such  that  if  t  >  0, 


(t)  at  the  time  t. 
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'f] ^  i3  the  vector  function  such  that  if  t  0. 

V^it)  =  (Rjj(t)  •  (t))  (t). 

is  the  simple  graph  such  that  if  t  0, 

V^it)  =  v^it)  ijj  (t)  at  the  time  t. 

r  is  the  vector  function  such  that  if  t  0,  r(t)  is  the  vector  from  the 
point  Q  to  the  jig  condition  location  of  the  airframe  reference  point  N  at 
the  time  t. 

d  is  the  simple  graph  such  that  if  t  >1  0> 
d(t)  =  r(t)  *  at  the  time  t. 

s  is  the  simple  graph  such  that  if  t  0, 
s(t)  =  r(t)  *  I2  at  the  time  t. 

r  is  the  vector  function  such  that  if  t  >  0,  r„(t)  is  the  vector  from  the 
n  —  n 

point  Q  to  the  jig  condition  location  of  the  point  H  at  the  time  t. 

I„  is  the  vector  function  such  that  if  t  >  0,  Iu(t)  is  the  vector  from  the 

jig  condition  location  of  the  point  N  to  the  jig  condition  location  of  the 

point  H  at  the  time  t  and  each  of  1  b  =  1,  3  is  a  number  such  that 

H 

ij^(t)  =  ij^^’  ig  (t). 
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Fj^p  is  the  vector  function  such  that  if  t  0,  r^pCt)  is  the  vector  from  the 
point  Q  to  the  ooirt  HP  at  the  time  t. 

v„  is  the  vector  function  such  that  if  t  >  0, 
n  “ 

Vjj(t)  =  (ry(t)  *  f^) 

Vu  is  the  simple  graph  such  that  if  t  >  0, 
n  “ 

VnCt)  =  Vjj(t)  at  the  time  t. 

Ij^p  is  the  vector  function  such  that  if  t  0,  is  the  vector,  with 

magnitude  the  number  l^^p,  from  the  point  LP  to  the  point  HP  at  the  time  t. 

Ij^p  is  the  vector  function  such  that  if  t  0,  with 

magnitude  the  number  l^p,  from  the  jig  condition  location  of  the  point  H  to 
the  point  LP  at  the  time  t. 

is  the  vector  function  such  that  if  t  0,  ijjy(t)  is  the  vector,  with 
magnitude  the  number  1^^^,  from  the  point  H  to  the  point  HP  at  the  time  t. 

djjp  is  the  vector  function  such  that  if  t  0, 

dHp(t)  =  (rjjp(t)  i^)  i^. 

dj^p  is  the  simple  graph  such  that  if  t  2^  0, 
dj^p(t)  =  dj^p(t)  at  the  time  t. 
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v„„  is  the  vector  function  such  that  if  t  >  0, 

HP 

VHp(t)  =  ^3^  -^3- 

v„„  is  the  simple  graph  such  that  if  t  >  0, 

HP 

VHp(t)  =  Vjjp(t)  at  the  time  t. 

I„  ,  I„  and  iu  are  the  unit  vector  functions  such  that  if  t  >  0, 
I„  (t)  has  the  direction  of  -  I„„(t),  !„  (t)  is  i„  (t),  and 

H^  HP  ^2  ^2 

i„  (t)  is  orthogonal  to  I„  (t)  and  !„  (t)  and  is  directed  dovm 
^3  1  2 

with  respect  to  the  pilot  at  the  time  t. 


Fjjp  and  Fj^  are  the  vector  functions  and  F^^p  and  F^^  are  the  simple  graphs 
such  that  if  t  >  0,  the  force  on  the  point  HP  from  the  arresting  hook  is 


FHp(t)  = 

=  Fjjp(t)  (t)  =  -  Fj^(t)  (t)  at  the  time  t. 


is  the  vector  function  and  F„  and  F„  are  the  simple  graphs  such  that 

Hi  Hi  Hi 

if  t  >  0,  F„  is  the  component  of  the  vector  F„  such  that 
—  H  ^  H 

F„  (t)  =  (F„  (t)  *  i„  (t))  I„  (t) 

Hi  Hi  Hi  Hi 

A 

=  F  (t)  I„  (t)  =  F„  (R0(t))  i„  (t)  at  the  time  t. 

Hi  Hi  Hi  Hi 


Each  of  Fji  ;  a  =  1 ,  3  is  the  simple  graph  such  that  if  t  0, 
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F„(t)  =  F„^(t)  i  at  the  time  t. 
n  n  3 

T  is  tne  vector  function  such  that  if  c  G,  T,  is  tne  force  on  the 
point  HP  from  the  cable  segment  between  the  points  HP  and  SL  at  the  time  t. 

f_  is  the  vector  function  such  that  if  t  >  0,  T„(t)  is  the  force  on  the 

n  -“A 

point  HP  from  the  cable  segment  between  the  points  HP  and  SR  at  the  time  t. 

In  addition,  the  following  simple  graphs  are  needed  in  the  arresting 
force  equations: 

k_.„  is  the  simple  graph  such  that  if  t  >  0, 
oAn 

kBAp(t)  =  |lL(t)|  -  llp(t)|  at  the  time  t. 

is  the  simple  graph  such  that  if  t  is  in  [0,  t^], 
between  the  arresting  hook  shank  and  the  unit  vector  I  (t)  at  the  time  t. 

A  A 

F^  is  the  simple  graph  such  that  if  x  is  in  [0,  F^(x)  is  the  magnitude 

of  the  arresting  force  at  the  cable  runout  x. 

Q.  ^  is  the  simple  graph  such  that  if  t  0,  Q.^(t)  is  the  angle  measured 
clockwise  from  the  straight  line  that  includes  the  points  HP  and  SL  to  the 
straight  line  that  includes  the  points  SL  and  S  at  the  time  t. 

Q P  is  the  simple  graph  such  that  if  t  >  0,  is  the  angle  measured 

counterclockwise  from  the  straight  line  that  includes  the  points  HP  and  SR 


is  the  angle 
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to  the  straight  line  that  includes  the  points  SR  and  S  at  the  time  t. 

is  the  simple  graph  such  that  if  t  0,  2  €^(t)  is  the  angle  measured 
counterclockwise  from  the  straight  line  that  includes  the  points  HP  and  SL 
to  the  straight  line  that  includes  the  points  HP  and  SR  at  the  time  t. 

The  cable  bisector  line  is  the  straight  line  that  is  in  the  plane  of  the 
cable  and  contains  the  point  HP  and  at  the  time  t,  is  the  angle 

measured  clockwise  from  this  line  to  the  straight  line  containing  the  points 
SL  and  HP. 

^  ^  is  the  simple  graph  such  that  if  t  0,  is  the  angle  measured 

clockwise  from  the  cable  bisector  line  to  the  straight  line  that  includes 
the  points  LP  and  HP  at  the  time  t. 

is  the  simple  graph  such  that  if  t  >  0,  €u  (t)  is  the  largest  number 

for  the  magnitude  of  the  angle  €  (t)  for  which  the  hook  will  not  slip  on 

n 

the  cable  at  the  time  t. 

A  A 

is  the  simple  graph  such  that  if  t  0,  has  the  magnitude 

1 6  (t)|  and  the  sign  of  €  (t)  at  the  time  t. 

Hl  H 

Q  ^  is  the  simple  graph  such  that  if  t  ^  0,  Cl^(t)  is  the  angle,  measured 
positive  from  the  ground  plane,  that  is  formed  by  the  intersection  of  the 
ground  plane,  the  cable  plane  and  plane  which  is  normal  to  both  the  ground 
and  cable  planes  at  the  time  t. 
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If  the  effects  of  the  airframe  deformation  are  neglected,  the 


orientation  of  the  hook  at  the  time  t  a  0  may  be  represented  by  the  diagram 
below.  In  this  diagram  the  point  N,  H  and  HP  are  coplanar. 


N 


Therefore,  it  is  seen  that  if 

I„(0)  =  1„»  r  (0), 

D 

then 

ll„p(0)  =  ?„p(0)  •  I,  =  d(0)  ♦  a/  -  S'**  IfiV  ' 

.  8^”  l„v  >> 

0  b 

If  the  hook  is  above  the  ground  for  t  <  t^,  then  the  diagram  below  may 
be  used  to  orient  the  hook. 
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I 


3 


ground  surface 


HC 


From  the  diagram  above  it  is  seen  that 

rjjp(t)  =  rj^(t)  + 

and  it  follows  that 

djjpCt)  =  rjjpCt)  ■  =  d(t)  +  dp'’  -  5"’  Ipp  sin{(f>^  ) 

*  8^"  l„v 

0  b 


If  the  hook  is  in  contact  with  the  ground  at  a  time  t  <  t  then  the 

angle  (p>^(t)  is  determined  as  follows:  Suppose  that  each  of  a  and  b  is  a 

simple  graph  such  that  if  t  is  in  [0,  t„],  the  vector  I„,,(t)  may  be 

C/  H  V 

expressed  by 


luu(t)  =  a(t)  i  (t)  +  b(t)  r  (t)  at  the  time  t, 
HV 
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and  if  it  assumed  that  the  point  HP  remains  a  fixed  distance  above  the 
ground  until  the  cable  is  picked  up  by  the  arresting  hook,  that  is 


VRpCt)  - 


VRCt)  ^3  “  ~  ^3  “  ''hP^‘^C^  ^3  "  Tb, 

b 


Thus,  it  is  evident  that 


VRCt)  =  -  a(t)  Tg^Ct)  -  b(t)  -  Vjjp(t^), 


1  (t;|^  »  a(t)^  +  b(t)^  = 


Consequently, 


a'(t)  [(7yT~^  ^  t— ^Vl - T~  ^  - 2 - ’ 

VRCt)  2  Vj^(t)  Vj^p(t^) 

u  B  (  (5  B 


-  (v„(t)  +  v^p(t|,))  -  a(t)  "yj’ct) 

77^^  ^ 


Therefore , 
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tan(C^^(t)) 


a(t) 


b(t) 

The  position  of  the  hook  at  the  time  of  cable  engagement  with  the  hook 
(i.e.  at  time  t  =  t^,)  is  shown  below. 


At  this  time  the  point  HP  is  at  the  height  v  (t  )  above  the  ground  and 

nr 

also 


“^BAR 


The  number  d_._  is  defined  as 
BAR 

dnpCtc)  - 


After  the  cable  has  made  engagement  with  the  hook,  it  may  now  pivot 
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abO'-.t  the  lateral  pivot  point  LP.  The  simple  graph  may  now  be 
determined  from  the  geometric  arrangement  shown  in  the  diagram  below 


It  is  seen  that 

(iLp(t)  »  S„(t)  •  I3  =  -  .„(t)  -  .3, 

and 

dhpCt)  +  i-|  =  "  ^1  ~ 


Also,  it  is  supposed  that 


(Ig  (t)  X  ^ 


Thus,  with 
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=  ^LP’ 

and 

dj^(t)  =  r^it)  •  i^  =  d(t)  + 

b 

I^p(t)  =  a^(t)  Ig  (t)  +  b^(t)  Ig  (t), 
^.,(t)  =  a^Ct)  i^(t)  +  bpCt)  iy 
it  follows  that 


a^(t)  +  b^(t)  +  b^Ct)  =  -  Vg(t)  -  Vg,  (F-1) 

a^(t)  TgVt)  +  b^(t)  ■"  ^2^^^  =  ‘ 

a^^(t)  +  b^^(t)  =  l^p^,  (F-3) 

-  ^gVt)  a^(t)  apCt)  +  Ifg^t)  b^(t)  ap(t)  (F-4) 


-  ^g^(t)  a^(t)  bpCt)  +  yg^(t)  b^(t)  bp(t)  =  0, 

=  ap^(t)  +  bp^Ct).  (F-5) 

The  following  procedure  is  used  for  solution  of  these  equations: 
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(1)  Assume  a  number  for  b^(t). 

(2)  Calculate  a^(t)  from  equation  (F-3)  (it  is  assumed  that  a^(t)  is 
negative) . 

(3)  Calculate  b^Ct)  from  equation  (F-1). 

(4)  Calculate  a2(t)  from  equation  (F-2). 

(5)  Calculate  a  new  candidate  for  b^(t)  from  equation  (F-4). 

(6)  Go  to  step  (2)  of  this  process. 

This  procedure  is  terminated  when  the  change  in  b^(t)  is  within  the 
required  precision. 

The  simple  graph  T/  is  determined  from  the  equation 

*  n 

V^it)  =  (?(t)  +  ij^(t)  +  I^pCt))  • 

from  which  the  following  expression  is  derived: 

=  s(t)  -t-  Ij^'^Ct)  7/(t)  +  a^(t)  ^B^Ct)  +  b^(t) 

b  1  3 

The  cable  geometry  at  a  time  t  in  the  plane  of  the  I„  (t)  and  I„  (t) 
vectors  is  shown  in  the  diagram  below. 
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Therefore,  since 

<V'>  * 

(I„p(i)  *  I,(t))  ■  =  -  S„ct), 

(I^p(t)  +  I^(t))  •  i^^(t)  =  -  (be  +  '^pjCt)), 

(I^p(t)  +  IpCt))  •  (t)  =  (be  -  ■^y(t)), 

it  follows  that  after  the  indicated  vector  operations  are  performed, 

l^p  3in(Q^(t)  +  ^^(t)  -  ^fj(t))  +  sin(a^(t))  =  ^pj(t),  (F-6) 

l^p  sin(a^^(t)  +  €^(t)  -  sin(ap(t))  =  ^^(t),  (F-7) 
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Ijjp  cosCQj^Ct)  +  €^(t)  -  €jj(t))  ♦  oo«(aj^(t))  =  be  ♦  (F-8) 

-  l^p  cosCQj^Ct)  +  €^(t)  -  €jj(t))  +  lp(t)  oo«(ap(t))  =  be  -  (F-9) 

From  equations  (F-6)  and  (F-8)  it  la  seen  that 
2  (bo  +  l^Ct)  oos(a^^(t))  +  2  ^^(t)  Ij^(t)  8in(a^^(t)) 

+  (l^p^  -  ^/(t)  -  (be  +  '  0* 

and  with 

x(t)  =  2  (bo  +  ^jj(t)) 
y(t)  =  2  ^jj(t) 

z(t)  =  (Ijjp^  -  ^/(t)  -  (bo  + 
it  follows  that 

x(t)  oos(Q^(t))  +  y(t)  sin(Q^(t))  x(t)  *  0, 
and  oonsequently 

(x^(t)  +  y^(t))  «in^(a^^(t))  +  2  y(t)  z(t)  sin(aL(t))  (F-10) 

+  (z^(t)  -  x^(t))  =  0. 

If  equation  (F-6)  la  substituted  into  equation  (F»7)  and  equation  (F-8) 
is  substituted  into  equation  (F-9)  it  ia  found  that 

l^(t)  ain(a^(t))  -  Ij^(t)  sln(Qp(t))  m  0, 

l^(t)  003(Q^^(t))  +  lp(t)  oos(ap(t))  =  2  bo. 
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Thus , 


tan(Clj^(t) ) 


sin(Q^(t) ) 


— -  C03(CL  (t) ) 
L(t)  ^ 


and 


(F-11) 


l^(t)  sin(Q^(t)) 

l„(t)  =  - 

sin(Q  j^(t) ) 


(F-12) 


Suppose  that  |6„  (t)!  (the  limit  hook  angle  with  respect  to  the  cable 
bisector  at  the  time  t)  is  the  largest  number  for  |€^(t)j  for  which  the 
hook  will  not  slip  on  the  cable. 


Therefore,  if 

l€^(t)l  is  less  than  or  equal  to  (t)|, 

L* 

then  the  hook  will  not  slip  on  the  cable  and  consequently 

-  l„(t)  .  k3,„(t).  (F-.3) 

The  numbers  1,  (t),  l^Ct),  Q,(t)  and  Cln(t)  for  this  case  may  be 
L  n  L  n 

obtained  as  follows: 

(1)  Assume  a  number  for  l^(t). 

(2)  Calculate  Q^^Ct)  from  equation  (F-10). 

(3)  Calculate  Q„(t)  from  equation  (F-ll). 

n 
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(4)  Calculate  l_(t)  from  equation  (F-12). 

n 

(5)  Calculate  a  new  candidate  for  l^(t)  from  equation  (F-13). 

(6)  Go  to  step  (2)  in  this  process. 


This  procedure  is  terminated  when  the  change  in  is  within  the 

desired  precision. 


The  angle  €„(t)  may  be  determined  from  the  equation 


6 


a^it) 

2 


a,(t) 

2 


(F-14) 


Further,  the  angle  €„(t)  may  be  found  by  combining  equations  (F-7)  and 

n 

(F-8)  above  as  follows: 


Equations  (F-6)  and  (F-8)  may  be  rearranged  to  determine 

^jj(t)  -  l^^t)  sin(Cl^(t)) 

sin(Q^(t)  +  ^^(t)  -  '  - » 

^HP 

and 

be  +  cos(Cl^(t)) 

cos(Q^(t)  +  €^(t)  -  €jj(t))  =  - 

^HP 


Therefore, 

sin(€j^(t))  =  sin(Q^(t)  +  €^(t))  cos(Q^(t)  +  -  ^^^(t)) 

-  cos(aj^(t)  +  sin(a^^(t)  + 
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co3(€jj(t))  =  cos(Qj^(t)  +  6(,(t))  cos(aj^(t)  +  €^(t)  - 
+  sin(Q^(t)  +  €^(t))  sin(a^(t)  +  ^^(t)  -  6jj(t)), 
and 

sin( ^  yCt) 

tan(€  (t))  =  - . 

cos(  €  jj(t) ) 

In  the  event  that  |€j^(t)|  as  found  from  the  equation  above  is  greater 
than  the  limit  angle  |€„  (t)j  then  the  following  definition  is  needed: 

A 

Suppose  that  €  is  a  simple  graph  such  that  if  t  >  0, 

H  “ 

€„(t) 

A  H 

€  (t)  =  - ,  €„  (t)  at  the  time  t. 

“  |€„(t)l  "l 

A 

With  €  (t)  substituted  for  €  (t),  equations  (F-6)  and  (F-8)  above  may 
n  n 

be  combined  to  derive  the  following  equation  for  the  number  1  (t): 

+  (be  +  cos(Q^(t)  +  ^^Ct)  -  €j^(t)))  . 

The  simple  graph  1^  may  be  determined  from  the  following  procedure: 

(1)  Assume  a  number  for  1  (t). 

Li 

(2)  Calculate  ([l^(t)  from  equation  (F-10). 

(3)  Calculate  Q  (t)  from  equation  (F-11). 

n 

(4)  Calculate  l„(t)  from  equation  (F-12). 

n 
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(5)  Calculate  €  (t)  from  equation  (F-14). 

u 

(6)  Calculate  Ij^(t)  from  equation  (F-15). 

(7)  Go  to  step  (2)  in  this  process. 

This  procedure  is  terminated  when  the  change  in  Ij^(t)  is  within  the 
desired  precision. 

The  cable  tensions  and  hook  forces  may  be  computed  from  an  examination 
of  the  forces  on  the  point  HP  as  shown  in  the  diagram  below. 


From  this  Hagram  it  is  found  that 
FHp(t)  +  f^^(t)  +  fj^(t)  =  0, 
or 

|Fj^p(t)|  cos(€^(t))  =  |Tp(t)|  cos(  €^(t)  +  lf^^(t)|  cos(€^(t)), 

-  I^HpCt)!  sin(€j^(t))  =  |fp(t)|  sin(€y(t))  -  |fj^(t)|  sin(€(,(t)). 
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With  the  friction  coefficient  between  the  hook  and  the  cable  the 
limiting  condition  for  the  hook  to  slip  on  the  cable  is  determined  from 

|f^(t)|  =  |fp(t)|  exp(/J.^  (  TT-  2 

=  |fp(t)|  exp(jLL^  CLp(t))) 

for  the  case  when  €„(t)  is  greater  than  zero  and  from 

H 

|fp(t)|  =  |f^^(t)l  expifl^  (  TT-  2  €^(t))) 

=  lT^(t)l  expifl^ 

for  the  case  when  ^  ^^(t)  is  less  than  zero. 


It  is  supposed  that  at  a  time  t  greater  than  t^,  there  is  a  component 

of  the  vector  F„  defined  as  F„  such  that 
H  n  ^ 

F  (t)  =  (t)  F  (t), 

Hi  Hi  Hi 


and  it  is  further  supposed  that  with 


RO(t)  =  lp(t)  sin(CLp(t)), 

A 

there  is  a  prescribed  simple  graph  F„  such  that 

F„  (RO(t))  =  F„  (t). 


From  this  definition  it  follows  that 

=  lTp(t)|  sinCQpCt))  +  iT^Ct)]  sinCa^Ct)). 


In  the  case  where  the  cable  is  slipping  on  the  hook  and  ^„(t)  is 

H 
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greater  than  zero, 


F  (t) 

_ ^ _ 

(exp(/Xc  +  QpCt)))  3in(aj^(t)  +  sinCQj^Ct))’ 


For  the  case  where  the  cable  is  slipping  on  the  hook  and  €„(t)  is  less 

H 

than  zero, 

F  (t) 

I  Vt)l  -  - ^ - . 

(exp(|LLc  +  QpCt)))  sin(ap(t))  +  sinCQ^^Ct))) 

Thus,  the  magnitude  of  the  axial  force  in  the  arresting  hook  shank  is 
=  !TR(t)|2  +  |TL(t)p  +  2  |fj^(t)|  |fj^(t)|  cos(?  €^(t)), 

and  the  number  €„  (t)  is  calculated  from 

“l 

(lf^(t)|  -  |fj,(t)|)  sin(  €j,(t)) 

sin(  €  (t))  = - . 

^  If„p<t)| 


In  the  case  where  the  cable  is  not  slipping  on  the  hook. 


iFHp(t)  =  Fj^p(t) 


F  (t) 


HP' 


sin(a^(t)  + 


|F^p(t)|  sin(  €jj(t)  + 
sin(2  ^(,(t)) 


158 


Fh  (t)  -  |T^(t)|  sin(a^(t)) 

It  (t)|  =  — - . 

^  sinCQpCt)) 

The  diagram  below  shows  the  relationship  at  the  time  t  between  the  unit 

vectors  I„  (t);  b  =  1,  3  and  the  unit  vectors  I  ;  c  =  1,  3. 

H,  0 


From  this  diagram  it  is  seen  that 
(t)  =  cos(Q^(t))  -  sin(Q^(t)) 

(t)  =  sin(a^(t))  +  cos(Q|^(t)) 

where 

-  b^Ct) 

tan(Qp(t))  =  - . 

a^Ct) 

This  transformation  may  be  written  in  the  form 

~  ^b’ 

a  a 


59 


where 


=  C03(Q^(t)), 

=  -  sinCQ^Ct)) , 

Tjj^Ct)  =  1, 

2 

y^\t)  =  3in(a(,(t)), 

=  co3(a^(t)) , 
and  all  other  Q^{t)  -  0. 

The  tranaformauion  from  the  the  unit  vectors  f  (t);  a  =  1,  3  to  the 

a 

unit  vectors  I„  (t);  b  =  1,  3  may  be  expressed  by 
^b 

Ih  (t)  =  apj^(t)  (t). 
b  "b  "a 

From  the  definition  of  the  unit  vector  function  I  it  follows  that 
(t)  =  sin(a^^(t)  +  €j,(t)  -  ^^Ct))  Ijj  (t) 

+  cosCa^^Ct)  + 
where 

Qj^Vt)  =  sin(a^(t)  +  €^(t)  -  €j^(t)), 

=  cos(aj^(t)  +  €^(t)  -  €j^(t)), 

a„^(t)  =  0. 
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Therefore , 


i  (t)  =a„V)  T/(t)  i^, 

b  b  a 


The  unit  vector  function  i„  may  be  determined  from  the  vector 

”3 


equations 


(t)  •  I  (t)  =  0, 

®2  ”3 


iLp(t)  •  I„^(t)  =  0. 


(F-16) 


(a/t)  (t)  +  b,(t)  (t))  •  I„  (t), 

1  1 


(F-17) 


i„  (t)  •  I„  (t)  =  1. 

H3  H3 


(F-18) 


After  the  vector  operations  in  equations  (F-I6),  (F-17)  and  (F-I8)  have 
been  performed,  they  may  be  rewritten  as  follows: 


23d 

(a^(t)  +  b^(t)  Tg^'Ct))  Tg^t)  = 

1  23d 

a/wa„\t)  5,,  =  1. 


These  equations  may  be  used  to  find  the  direction  cosines 

Qu^Ct)  ;  a  =  1  ,  3. 

3 


From  the  definition  of  the  vector  function  F  , 

H 
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■  “  ^HP' 

In  addition,  F„  may  be  determined  from 
n 

FnCt)  =  (t) 

where 

F^\t)  =  -  f^pCt)  Gos(a^(t))  sin(Q^(t)  +  ^^Ct)  -  ^jj(t)), 

F^^(t)  =  -  f'jjp(b)  cosCQ^Ct)  + 

Fj^^(t)  =  F^pCt)  sin(Q^(t))  3in(Q.^(t)  +  ^  ^it)  -  €j^(t)). 

Since  these  relations  are  true  for  each  number  t  greater  than  t^,  it 
follows  that 

Fh^  =  -  F^p  cosLOL^]  sintQj^  +  -  €y], 

V  =  ■  ''hP  -  ^h'- 

'  ^HP  3in[Q  ^ 
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APPENDIX  G 


TRANSFORMATIONS 


The  orientation  of  the  airframe  relative  to  the  ground  is  shown  in  the 
diagram  below.  The  point  Q  is  the  ground  reference  point  and  the  point  N  is 
the  airframe  reference  point. 


The  transformation  from  the  ground  based  unit  vectors  (I^;  b  =  1,  3)  to 

the  airframe  fixed  unit  vectors  (1^  ;  c  =  1,  3)  may  be  established  through 

c 

the  use  of  the  yaw,  pitch  and  roll  Euler  angles.  This  may  be  accomplished 
with  the  aid  of  the  diagram  below. 
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To  illustrate  how  these  Euler  angles  are  used  in  the  transformation 
from  the  ground  based  unit  vectors  to  the  airframe  fixed  unit  vectors 


suppose  that  if  t  >  0,  the  airframe  unit  vectors  i  (t);  b  =  1,  3  are 


initially  oriented  the  same  as  the  ground  based  unit  vectors  c  =  1,  3. 


Now  suppose  that  the  unit  vectors  i  (t)  and  (t)  are  rotated  as  shown  by 
/  ^  _  2 

the  "yaw  angle"  Y {t) ,  with  the  unit  vector  (t)  held  fixed.  In  this 

®3 


position  it  is  seen  that 


i„  (t)  =  cos( 


V^(t))  +  sin(V^(t)) 


i  (t)  =  -  BiniYit))  I  +  cos( 


ig^(t)  = 
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From  the  orientation  described  above  now  suppose  that  the  airframe  unit 
vectors  (t)  and  (t)  are  rotated  as  shown  by  the  "pitch  angle"  Q  {t) , 

with  the  unit  vector  (t)  held  fixed.  In  this  position  it  is  found  that 


i„  (t)  =  cos 


(6^(t))  cos(^(t))  +  cos(^(t))  sin(*l^(t))  -  sin(C  (t)) 


(t)  =  -  sin('^(t))  i,  +  cos(^(t))  i-. 

82  ^  2 

ig  (t)  =  sin((?(t))  cos(^(t))  +  sin(6*(t))  sin(V^(t))  +  cos((9(t))  I, 


Finally,  suppose  that  the  airframe  unit  vectors  I_  (t)  and  (t) 

®2 

are  rotated  as  shown  by  the  "roll  angle"  (pit),  with  the  vector  (t) 

^  8^ 

held  fixed.  In  this  position  the  airframe  unit  vectors  are  aligned  with  the 
airframe  as  defined  in  the  main  body  of  the  report  and  the  transformation 
from  the  ground  based  unit  vectors  to  the  airframe  unit  vectors  is  found  to 


ip  (t)  =  cos(!7(t))  cos(  ^(t))  I,  +  cos  iOit))  sin  (V^(t)) 

Bi  ’  2 

-  sin(  (9(t) ) 

Ig  (t)  =  (-  cos(c/:(t))  sin(^(t))  +  sin«^'(t))  sin((5(t))  cos(V^(t))) 
+  (cos(<'/j  (t) )  cos(^(t))  +  sin«^(t))  sin(^(t))  sin(V^(t))) 

+  sin((^^(t))  cos((9(t))  1^, 

Ig  (t)  =  (sinCi^Kt))  sin(V^(t))  +  cos((^(t))  sin  (t^(  t))  cos(  Pit))) 
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+  (-  sin((^(t))  cos(  ^(t))  +  co»(C^(t))  8in(6^(t))  sin(^(t))) 

+  cos(cf  (t))  cos(6’(t)) 

p  h 

Now  suppose  that  i  „  ;  c  =  1 ,  3  i»  a  simple  surface  such  that 

D 

c 

the  transformation  above  may  be  written  in  the  form 

ig  (t)  =  ^](t)  (G-1) 

c  c 

or 

‘b  "•>  =  <“-2) 

c  c 

In  the  Kinetic  Energy  Formulation  r  was  defined  as  the  vector  function 
such  that  if  t  0,  r(t)  is  the  vector  from  the  point  Q  to  the  jig  condition 
location  of  the  point  N  at  the  time  t.  Further,  r(t)  was  defined  as 

?(t)  =  ^^(t)  1^. 

It  follows  then  that 

r’(t)  =  ^‘’’(t)  =  Vg^’ct)  Ip  (t),  (G-3) 

b 

where  V  ^(t)  was  defined  as  the  bth  body  axis  component  of  the  velocity 

D 

vector. 

From  equation  (G-3)  it  is  seen  that 

C'*'  Sbd  =  ■  ^B 

d 


166 


But,  from  equation  (G-2)  it  is  found  that 


(I,  •  Ib  <«>  = 

d  d 

Therefore, 

VB^t)  =  X^t)  S“  ^  = 
d 

In  tHo  ?f;inet^.lc  Energy  Formulation  the  transformation  from  the 
generalized  coordinate  velocities  to  the  quasi-coordinate  velocities  was 
given  as 


It  follows  therefore  for  a  in  the  interval  [1,  33  and  b  in  the  interval 
[1,  3]  that 

=  V 

c 

The  transformation  from  the  quasi-coordinate  velocities  to  the 
generalized  coordinate  velocities  was  given  by 


Therefore, 


/S 


a 


b 


for  a  in  [ 1 ,  33,  b  in  [ 1 , 


33  and  c  in  [1,  33  it  follows  that 
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The  terms  where  a  ie  in  [4,  6]  and  b  is  in  [4,  6]  may  be  determined 

as  follows:  Since  the  relation  in  (G-1)  is  true  for  each  positive  number  t 

then  the  vector  function  Ig  may  be  differentiated  to  obtain 

c 


I.;  •  »■  ■ 


But , 


SfhS'”'  ■ 

d 


and 


'  ■  '^B  *  ^B 

c  c 


=  e 


bdf 


«b‘ 


Sof  Sd. 


Consequently, 


bd 


(G-4) 


When  the  operations  in  equation  (G-4)  are  performed  it  is  found  that 

=  (^  '  -  sinC^’]  ^  ' . 

Q/  =  coM[(p]  0  '  +  cos  [01  sin[c/->] 

=  -  sin[C^]  0'  +  cos [6']  cos 
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In  many  cases  it  is  desirable  to  have  a  reference  system  that  is  fixed 

relative  to  the  airframe  but  rotated  to  be  aligned  with  a  gear  component 

such  as  a  shock  strut.  Such  a  reference  system  can  be  defined  through  the 

use  of  the  Euler  angles  Tj  ^  and  defined  in  the  main  body  of  the  report 

a  a 

and  illustrated  in  the  diagram  below. 
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170 


From  this  orientation  now  suppose  that  the  a  gear  unit  vectors  (t) 

ll 


and  i„  (t)  are  rotated  to  their  final  position  as  shown  in  the  diagram 
G 

r 

below  through  the  roll  Euler  angle  »  with  i„  (t)  held  fixed. 

U  U 


In  this  alignment  it  is  found  that 


i  (t)  =  cos('77_  )  Ig  -  siniTJ  )  I  (t). 


G  '  B, 
a  1 


G  '  B. 
a  3 


I_  (t)  =  sin 


Si  1 


+  cos(';:,Q  )  ig  (t) 
a  2 


+  cosiTj ^  )  sin(^Q  )  Ig  (t), 

a  3 


Ig  (t)  =  sin(77g  )  cos(  )  I„  (t)  -  sin(^^  )  (t) 


-a  ^ 


‘^a  "2 


+  cos(7^g  )  cosC^g  )  Ig  (t), 

^a  “3 


which  may  be  written  in  the  form 


Ig  --X'’  h 

a  a  b 

c  c 


(G-5) 


The  i„  (t);  c  =  1,  3  unit  vectors  are  illustrated  in  the  diagram 

c 

a 

c 

below  for  the  case  where  it  is  desired  to  have  the  shock  strut  piston 


centerline  aligned  a  time  t  with  the  vector  I_  (t). 

U 

^3 
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/ 

a 

For  gears  with  a  castoring  degree  of  freedom  it  is  convenient  to  have  a 

set  of  unit  vectors  fixed  in  the  gear.  These  unit  vectors, 

3q  (t);  b  =  1,  3  are  shown  relative  to  the  (t);  b  =  1,  3  unit  vectors 
%  ®b 

in  the  diagram  below. 


Jq  (t) 
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It  is  seen  that 


(t)  =  cos  (Or  (t))  Jr  (t)  -  sin(^^  (t))  (t). 


G  '  "  ■'G 
a  a. 


a  a. 


(t)  =  sin(^g  (t))  (t)  +  cos(^g  (t))  (t), 

a„  a  a,  a  a„ 


ic  (t)  =  (t). 


which  may  be  expressed  as 


ig  (t)  =  Qg  ^(t)  (t), 

a  a  a, 

c  c  b 


Now  suppose  that  if  t  0,  c  is  in  [  1 ,  3]  and  d  is  in  [1, 

X, 


there  is  a  transformation  ''(t)  such  that 

Cj 

^d 

(t)  =  Xg  ‘^(t)  (t)  at  the  time  t. 


It  is  evident  that 


J'c  ■  ‘c  Scg-ac”'')  S,,. 

a  a  a  a 

d  g  d  g 


Thus  , 


Xc  '<'>  -  Qc  Sd£  S'®' 


Therefore,  it  is  found  that 


Jo  -^G  Sg,  S'®  ic  (t). 


(G-6) 


3]  that 


(G-7) 
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Equations  (G-5)  and  (G-7)  may  be  combined  to  obtain  the  transformation 


3o  Sde 

a .  a  a  b 

d  g  c 


Thus,  if 


a„  Nt)  "  X  ^t)  Sh„  S°®. 


G  t^G 
a  a 
g  G 


3q  (t)  =  Qq  °(t)  Ig  (t). 
^d  ^d 


(G-8) 


Equation  (G-8)  may  be  rewritten  in  an  expanded  form  as 

3q  (t)  =  (cos(^Q  (t))  gos(7^q  )  +  sin(;^Q  (t))  sinCT^^  )  sin(^Q  ))  Ig  (t) 
a^  a  a  a  ^a  a  “l 

♦  ain(^  (t))  co3(^„  )  I  (t) 


'^a  «2 


+  (-  cos  (0Q  (t))  sin  ^Vg  )  +  sin(^Q  (t))  co3(7^q  )  3in(^Q  ))  Ig  (t), 


a  3 


3q  (t)  =  (-  sinC^Q  (t))  cosCT^q  )  +  cos(^Q  (t))  sinCT^^  )  sini^^  ))  Ig  (t) 


a  1 


+  cos  (t))  cos  )  Ig  (t) 

a  a  2 


+  (sinC^Q  (t))  sinCT^Q  )  +  cosiO^  (t))  co 3 {Tj ^  )  3in(^^  ))  Ig  (t), 


a  3 


3q  (t)  =  3in(7^^  )  cos(^|,  )  Ig  (t)  -  3in(^  '  I  (t) 

a^  a  a  1  °2 

+  cos  (7^  Q  )  cos(  <«• 

a  ^a  3 


The  a  gear  wheel  fixed  unit  vectors  are  oriented  as  shown  in  the 
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diagram  below 


'  axle 
centerline 


At  a  time  t  the  vector  k-  (t)  is  aligned  with  the  axle  centerline 

U 

^2 

and  is  expressed  in  terms  of  the  a  gear  fixed  unit  vectors 


(t);  b  :  1,  3  by  the  relation 

U 

^b 

(L)  =  3^  ^(t)  (t), 


(G-9) 


where  the  (t);  b  =  1,  3  are  direction  cosines  determined  from  the 

U 

^2 

rotation  of  the  axle  that  is  due  to  stroking  of  the  a  gear  shock  strut.  The 
general  relationship  between  these  two  sets  of  unit  vectors  is 


(t)  =  ^(t)  L  (t), 


G  ' 

a  a 

c  ^b 


where  the  unit  vectors  (t)  and  (t)  are  orthogonal  to  each  other  and 

Si  ■’ao 
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to  k-  (t),  but  otherwise  arbitrary.  For  the  case  where  the  a  gear 
G 

^2 

castors,  it  is  usual  for  k_  (t)  to  be  orthogonal  to  (t). 

u  G 


The  wheel-ground  reference  system  unit  vector  functions 
;  b  =  1,  3  are  defined  in  Appendix  E  in  terms  of  the  vector  functions 

k^  and  b  =  1 ,  3- 
^2 

At  a  time  t  the  transformations  du^(t)  and  Q.,,°(t),  which  are  defined 

n  H  , 

c  d 

in  Appendix  F,  are  used  to  relate  the  unit  vectors  b  =  1,  3, 

T„  (t);  c  =  1,  3  and  I  (t);  d  =  1,  3  through  the  relations 
H  H  , 

c  d 


^b’ 
c 


I  (t)  =  a  ^^(t)  (t), 

”d  "d  ”c 


I  (t)  =  a  ®(t)  ^ 


H  *b' 

e 
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